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1.  INTRODUCTION  AND  SUMMARY 


This  report  summarizes  the  study  of  a stratified 
dielectric  slab  spatial  filtering  technique  to  provide 
element  pattern  control  on  large  planar  arrays.  Specific- 
ally, the  objectives  of  the  study  were  to: 

- Develop  a computer  optimization  technique  for  synthesis 
of  stratified  dielectric  slab  spatial  filters  subject 
to  realistic  constraints. 

- Synthesize  optimum  filter  designs  subject  to  required 
rejection  characteristics  and  frequency  bandwidth 
requirements  at  X-band. 

The  stratified  dielectric  slab  spauial  filtering 
technique  was  first  described  by  Maillouxv  , wherein 
the  properties  of  filtering  in  the  spatial  dcma: n were 
outlined  and  the  analytical  synthesis  of  filters  possessing 
Chebyshev  characteristics  was  presented.  The  present 
work  extends  Mailloux's  earlier  results  by  providing  a 
numerical  technique  for  synthesizing  filters  with  arbitrary 
spatial  characteristics.  As  a consequence  of  the  method 
developed  here,  the  physical  constraints  encountered  in 
practical  filter  design  (e.g.,  material  availability  and 
layer  thickness)  may  be  directly  parametrized,  thereby 
providing  a practical  filter  design  as  the  output  of  the 
synthesis . 

The  principle  point  of  departure  of  the  analytical 
approach  to  filter  synthesis,  as  typified  by  Mailloux's 
results,  and  the  numerical  approach  presented  here,  is 
associated  with  the  ultimate  design  goal.  In  the 
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analytic  approach,  the  degree  of  specification  is  limited, 
essentially,  by  the  algebraic  complexity  of  the  constituent 
equations,  with  the  result  that  the  design  goal  must 
necessarily  have  limited  scope  with  respect  to  the  overall 
degree  of  freedom  associated  with  the  physical  structure. 
For  example,  in  the  analytical  procedure  given  by  Mailloux, 
this  effective  restriction  results  in  the  specification 
of  the  filter  only  in  terms  of  the  desired  electrical 
performance.  To  be  sure,  this  may  lead  to  practical 
spatial  filter  designs.  However,  the  procedure  results 
in  a specification  of  dielectric  materials  which  may,  or 
may  not,  be  obtainable  in  practice.  In  effect,  no  material 
selection  control  is  provided  by  the  analytical  synthesis 
procedure  and  the  practicality  of  the  synthesized  filter 
depends  strongly  on  the  availability  of  the  materials 
so  determined. 

Using  the  numerical  optimization  technique,  the 
overall  design  goals  may  be  factored  into  the  procedure. 
With  respect  to  the  selection  of  dielectric  material, 
this  procedure  is  also  deterministic.  However,  as  a 
particular  goal,  we  may  ask  the  optimizer  to  select  one 
or  all  of  the  materials  from  a predetermined  list.  By 
judiciously  constructing  this  list,  such  that  only  avail- 
able materials  are  represented,  the  optimizer  will  always 
synthesize  a design  which  is  readily  fabricated.  By 
interference,  such  a design  would  be  low  cost  due  to  the 
specific  selection  of  "off  the  shelf"  materials. 

In  principle,  the  advantage  of  the  numerical  opti- 
mization technique  is  that  the  available  solution  set  is 
not  artificially  constrained  to  a limited  space,  spanned 
only  by  those  physical  parameters  which  are  ammenable  to 
direct  inclusion  in  an  analytical  expression.  It  is 
clear,  however,  that  in  practice,  not  all  physical  (elec- 
trical and  mechanical)  properties  of  the  structure  can 
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be  modeled  to  provide  a reasonable  optimization  parameter. 
In  particular,  during  the  course  of  this  study,  the 
design  parameters  material  density,  tensile  strength, 
sheer  strength,  dielectric  strength,  and  machinability 
were  found  to  be  of  second  order  importance  for  typical 
applications.  For  low  loss  microwave  materials  in  the 
range  l<er<25,  these  properties  either  varied  little,  or 
were  not  determinable  from  available  data.  Consequently, 
the  scope  of  the  synthesis  procedure  has  been  limited  to 
the  consideration  of  material  selection,  layer  thickness, 
and  frequency  bandwidth. 

The  synthesis  procedure  developed  in  this  study 
proceds  from  a wave  transmission  analysis  of  scattering 
of  incident  plane  waves  by  infinite  plane  stratified 
dielectric  slabs.  An  equivalent  circuit  representation 
of  the  layered  media  is  constructed,  and  from  the  scatter- 
ing matrix  description  of  the  constituent  networks  rep- 
resenting junction  discontinuities  and  length  of  trans- 
mission line,  a wave  transmission  formalism  is  obtained. 

A single  wave  transmission  matrix  A,  relating  the  incident 
and  reflected  traveling  wave  voltages  at  the  filter  input 
to  those  at  the  output,  is  constructed.  Since  the  power 
transmitted  with  respect  to  one  volt  at  the  input  is 
simply  (1/A^ ) , we  operate  directly  on  this  quantity, 
and  optimize  the  filter  response  with  respect  to  the 
desired  rejection  characteristics  over  the  frequency 
bandwidth.  The  optimization  algorithm  used  is  a modifica- 
tion of  the  RAZOR  search  technique  due  to  Bandler  and 
(2) 

MacDonald.  The  synthesis  procedure  is  given  in 

Section  2. 
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In  section  3,  an  eleven  layer  stratified  dielectric 
slab  spatial  filter  design  is  presented.  The  design  goal 
was  to  achieve  less  than  .1  db  insertion  loss  in  the  spatial 
pass  band,  and  greater  than  10  db  rejection  in  the  spatial 
stop  band  over  a 4%  frequency  band.  Since  the  permittivity 
of  most  microwave  materials  is  constant  from  S band  through 
X band,  the  layer  thicknesses  are  given  normalized  to  center 
frequency  wavelength.  In  addition,  the  designs  synthesized 
are  such  that  the  spatial  extent  of  spurious  passbands 
produced  for  TM  incidence  by  the  Brewster  angle  phenomenon 
is  minimized. 

Three  appendices  are  included.  Appendix  A summarizes 
the  modified  RAZOR  algorithm  and  associated  computer 
programs.  In  Appendix  B,  a table  of  available  low  loss 
microwave  dielectrics  is  presented.  This  table  gives 
a summary  of  microwave  materials  currently  available  and 
their  properties  (where  known) . Appendix  C is  a compilation 
of  computer  programs  developed  during  the  study. 
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2.  ANALYSIS  AND  SYNTHESIS  OF  PLANE  STRATIFIED 
DIELECTRIC  SLAB  SPATIAL  FILTERS 


The  plane  stratified  dielectric  slab  spatial  filter 

shown  in  Figure  1 is  configured  in  such  a manner  as  to 

provide  a multiplicative  element  pattern  modification 

of  the  type  shown  in  Figure  2 for  a large  planar  phased 

array  antenna.  The  filter  is  assumed  flush  with  the  array 

face.  Both  relative  permittivity  and  thickness  of  the 

slabs  are  arbitrary,  and  assumed  uniform.  All  media  are 

- 7 

assumed  to  have  permeability  uo(=4ttx10  henries/m). 

Basic  to  the  problem  of  synthesis  of  stratified 
dielectric  slab  spatial  filters  is  the  analysis  of 
arbitrary  transversely  infinite  layered  dielectric  sheets 
subject  to  plane  wave  incidence  with  TE  and  TM  polariza- 
tions. This  analysis  and  the  fundamental  limitations  of 
the  dielectric  slab  spatial  filter  are  presented  in  Section 
2.1.  The  synthesis  technique  is  presented  in  Section  2.2. 

A time  dependence  e-^  is  assumed  throughout. 

2 . 1 Analysis  of  Plane  Wave  Scattering  by  Plane  Stratified 
Dielectric  Sheets 

It  is  well  known  that  the  electromagnetic  fields 
in  the  vicinity  of  a large  planar  array  of  radiating 
apertures  may  be  decomposed  into  linearly  polarized  plane 
wave  constituents.  Consequently,  for  spatial  filter 
analysis  and  synthesis,  the  array  may  be  removed,  and  the 
properties  of  the  filter  may  be  examined  with  respect  to 
monochromatic  plane  wave  incidence  on  a transversely 
infinite  structure.  Since  arbitrary  linearly  polarized 
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Figure  1 - Dielectric  Slab  Spatial  Filter  over  Planar 

Phased  Array 
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sin  8 


Figure  2 - Multiplicative  Element  Pattern  Modifier  for 

Spatial  Filtering 

waves  in  the  vicinity  of  a dielectric  interface  may  be 
decomposed  into  TE  and  TM  components  with  respect  to  the 
surface  normal  which  remain  decoupled  upon  scattering, 
the  analysis  may  be  further  simplified  by  independently 
determining  the  transmission  and  reflection  character- 
istics of  the  device  for  the  individual  components.  The 
complete  solution  is  then  obtained  by  applying  super- 
position to  reconstitute  the  fields  in  all  regions. 

Consider  then  the  plane  stratified  dielectric 
sheet  shown  in  cross-section  in  Figure  3.  For  conven- 
ience and  without  loss  of  generality  plane  waves  are 
considered  incident  only  from  Z>0  at  an  angle  0 with 
respect  to  the  surface  outward  normal.  We  wish  to  deter- 
mine the  magnitude  and  phase  of  the  reflected  fields 
+ 

at  Z=0  , and  the  magnitude  of  the  transmitted  field  at 
Z=-D~  for  both  incident  polarizations  as  a function  of 
incidence  angle  0;  layer  permittivity , e 1 ; and  layer 
thickness,  t^. 
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INCIDENT  PLANE  WAVE 


Figure  3 - Stratified  Dielectric  Slab  Cross-Section 

Having  specialized  to  the  scattering  of  plane 
waves  which  remain  decoupled  at  the  dielectric  inter- 
faces, the  field  problem  may  now  be  reduced  to  an  equiv- 
alent transmission  line  problem  in  a standard  manner. 

The  equivalent  circuit  representation  is  shown  in  Figure 
4.  The  ith  dielectric  layer  of  the  physical  structure, 
with  relative  permittivity  er  and  thickness  t^ , is 

represented  by  the  i*"^1  transmission  line  segment  of  length 

t h 

t^.  The  characteristic  admittance  of  the  1 line  is  the 
wave  admittance  associated  with  the  dielectric  properties 
of  the  corresponding  layer  for  the  particular  incidence 
angle  6,  and  is  given  as 
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1) 


Y .= 
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(<i'u)Uo  for  TE  incidence 

we  . e /< . , for  TM  incidence 
no  1 


In  equation  (1),  w is  the  angular  frequency,  and  k . , the 

longitudinal  wave  number  in  the  l medium,  is  given  as 
0 


2) 


_W  / 2 

Ki  c sinQ 


since  the  filter  is  assumed  embedded  in  vacuum- c is  the 
speed  of  light  in  vacuum. 

It  is  a simple  matter  to  construct  the  voltage 
scattering  matrix  representation  of  the  circuit  in 
Figure  4.  Given  any  lossless,  reciprocal  two  port  net- 
work, as  shown  in  Figure  5,  the  voltage  scattering  matrix 
is  defined  by 


3) 


where  V and  V+  are  the  vector  representation  of  scattered 
and  incident  traveling  wave  voltages,  respectively.  By 
definition,  then,  the  elements  of  S are 


" Vi/Vl 
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Figure  4 - Equivalent  Circuit  for  Arbitrary  Layered 
Dielectric  Structure 


Figure  5 - Voltage  Labeling  at  Ports  of  Two  Port  Network 
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Each  transmission  line  segment,  and  each  junction  may  be 
represented  by  its  two  port  scattering  parameters.  By 
appropriate  manipulation  of  port  voltages,  a single 
scattering  matrix  may  be  obtained  for  the  network  equiv- 
alent of  the  circuit  in  Figure  4. 

A more  convenient  representation  of  the  circuit  is 
in  terms  of  its  traveling  wave  transmission  characteris- 
tics. The  traveling  wave  voltages  at  the  ports  of  a two 
port  network  are  related  by 


5) 


£ £2 


where  and  C_2  are  the  vector  representations  of  the 
incident  and  scattered  traveling  wave  voltages  at  ports  1 
and  2,  respectively,  as  defined  in  Figure  5.  By  properly 
associating  the  elements  of  and  C2  with  the  elements 
of  V and  V , the  elements  of  the  wave  transmission 
matrix,  A,  are  obtained  in  terms  of  the  elements  of  the 
voltage  scattering  matrix  S,  and  are  given  as 
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Representing  each  line  segment  and  junction  by  its 
two  port  wave  transmission  parameters,  Ad  and  A? , respec- 
tively, the  cascade  connection  of  networks  in  Figure  4 
results  in 


10)  C = A C- 

-o  — — fs 

where 


11) 


A = A] 


=o  id 

i=l 


is  the  complete  wave  transmission  representation  for  the 
network,  and  C^s  is  the  traveling  wave  voltage  vector 
at  Z = -D- . In  particular,  for  no  excitation  from  Z<-D 


12) 


Sfs 


and  the  transmission  coefficient  of  the  equivalent  net- 
work is  obtained  from  equation  (10)  as 


13) 


The  reflection  coefficient  at 


Z=0 


+ 


is 


then 


I 
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Formally,  equations  (13)  and  (14)  completely  specify 
the  scattering  properties  of  the  equivalent  circuit  in 
Figure  4 for  monochromatic  TE  or  TM  incidence,  and  it 
remains  only  to  determine  the  parameters  of  the  various 
junction  and  line  length  transmission  matrices.  These 
are  obtained  from  elementary  transmission  line  theory 
and  are  given  as 


15) 

A3  = A3  = (Y. 

111  122  1 

+vi+il/2Yi 

16) 

A3  = A3  = (Y. 

112  121  1 

-Yi+i»/2Yi 

for  the 

ifc^  junction,  and 

• 

17) 

ad  -i^i 

V = e 

18) 

4 = i/aJ 

22  X11 

19) 

„d  ^d  . 

Ai  = At  =0 

12  21 

for  the  i*"*"1  line  segment.  Examination  of  equation  (16) 
shows  that  for 


13 


20) 


Y _ Y 
i i + 1 


the  off-diagonal  matrix  elements  are  0.  Consequently, 
at  any  spatial  angle  for  which  equation  (20)  is  satis- 
fied, the  junction  is  transparent,  and  = 1.  From 

equations  (1)  and  (2),  this  equality  can  occur  only  for 
TM  incidence,  provided  ? ei+i" 

Equations  (1),  (2),  and  (20)  may  be  manipulated  to 

result  in  an  expression  for  the  incidence  angle,  0,  at 
which  a junction  is  transparent  to  TM  incidence.  This 
angle  is  given  as 


21) 


6T  = sin 


1+VEi.i 


where  f For  either  e^=l,  or  e^+^=l,  equation  (21) 

is  the  familiar  Brewster  angle  formula,  and  gives  the 

smallest  angle,  0=0,  for  which  the  junction  is  trans- 

r p 

parent.  For  both  / 1 and  + ^ ^ 1,  ©T  is  also  a 
spatial  angle  at  which  total  transmission  is  obtained 
at  the  junction,  provided. 

22)  ei+i  < 4/(ei-l) 


For  the  layered  dielectric  sheet  embedded  in  free  space, 
©T  is  bounded  by 
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23) 


sin”'W— <6  <90° 

V n 

where  c'  is  the  smaller  relative  permittivity  of  the 
first  and  last  dielectric  layers. 

As  a consequence  of  equation  (20),  the  layered 
dielectric  sheet  will  have  regions  of  high  transmission 
to  TM  incidence  at  angles  0>45°.  From  equation  (23), 
this  region  of  high  transmission  may  be  reduced  by 
selecting  the  dielectric  constants  of  the  outer  layers  to 
be  moderately  large.  A plot  of  the  left  hand  side  of 
equation  (23)  is  shown  in  Figure  6.  For  e'  = 3,  the 
lower  bound  in  equation  (23)  is  60°.  and  for  e'  = 5, 
the  lower  bound  is  65.9°.  For  e‘>5,  the  slope  of  the 
curve  is  very  slow,  and  little  advantage  is  gained  from 
the  larger  values. 


Figure  6 - Brewster's  Angle  versus  Relative  Permittivity,  e‘ 
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The  principle  task  in  stratified  dielectric  slab 
spatial  filter  synthesis  is  to  reduce  the  spatial  extent 
of  the  Brewster  angle  related  passbands  while  maximizing 
transmission  in  the  desired  spatial  passband  near  broad- 
side, as  shown  in  Figure  2. 

At  any  junction,  i,  of  the  equivalent  circuit,  the 
reflection  coefficient  is  given  as 


25) 


ri  ° Ai  /Ai 

12  X11 


= <Yi  -W/<Yi  +Vi> 


and  the  magnitude,  for  small  incidence  angles  0,  is  large 
for  permittivity  ratios  p =e ^/t  ^+^>>1  » or  P<<:1-  To 
achieve  spatial  filtering  it  is  therefore,  desireable 
to  select  geometries  which  give  permittivity  ratios  in 
these  ranges.  This  is  particularly  important  in  the  case 
of  TM  incidence  for  which  | T ? | is  a monotonically  de- 
creasing function  of  incidence  angle  out  to  6T.  Since 
the  range  of  available  microwave  dielectric  materials 
lies  in  the  approximate  limits,  l<e<25,  structures 
providing  the  spatial  filtering  characteristics  shown  in 
Figure  2 must  be  formed  by  alternating  layers  of  high 
>3)  and  low  (f^sr  1)  permittivity. 
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2 . 2 Synthesis  of  Plane  Stratified  Dielectric  Slab 
Spatial  Filters 

The  principle  objective  of  this  study  is  to  develop 
a technique  for  the  synthesis  of  practical  plane  strati- 
fied dielectric  slab  spatial  filters  for  use  with  large 
planar  arrays.  The  particular  technique  developed  is  one 
of  numerical  optimization  wherein  it  is  possible  to 
synthesize  filters  which  are  optimized  with  respect  to 
electrical  performance  in  both  the  spatial  and  frequency 
domains  while  constraining  the  solution  space  such  that 
only  practical  filter  configurations  result. 

The  spatial  filter  matrix  characterization  given 
in  equation  (11)  of  the  preceding  section  is  constructed 
such  that  all  electrical  and  physical  parameters  effect- 
ing filter  performance  are  modeled.  These  parameters  are 

- Plane  wave  polarization,  i.e.,  TE  or  TM  with 
respect  to  the  filter  normal 

- Plane  wave  incidence  angle,  6 

- Slab  relative  permittivity,  e^,  i=l,...,N 

- Slab  thickness,  t^ , i=l,...,N 

- Number  of  slabs,  N 

- Frequency,  f 


The  objective  of  the  synthesis  procedure  is  to 
determine  the  set  of  permittivities,  {e^},  and  slab 
thicknesses,  { t^  > , which  best  match  the  filter  perform- 
ance to  the  desired  performance  in  both  spatial  and 
frequency  domains. 
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The  desired  spatial  filter  performance  may  be 
specified  by  the  spatial  variation  of  the  power  trans- 
mission coefficient  as  shown  in  Figure  7.  In  the  figure, 
the  desired  performance  is  given  in  terms  of  forbidden 
regions  since,  in  the  spatial  domain,  specific  functional 
variation  control  of  the  measure  is  limited  by  physical 
realizability,  as  was  shown  in  Section  2.1.  In  the  region 
0<sin  0<sin  0^,  the  filter  is  to  be  maximally  transmissive, 
within  P^db  of  perfect  transmission.  In  the  region 
sin02<sin0<l , the  transmission  coefficient  is  to  be  less 
than  ?2db.  And  in  the  intervening  region,  sin0 , <sin0<sino 
the  transmission  is  to  be  monotonically  decreasing.  The 
forbidden  region  diagram  may  be  fixed  or  varying  in 
frequency . 


Figure  7 - Forbidden  Region  Diagram  of  Desired  Spatial 

Filter  Performance 
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For  incident  waves  which  are  TE  with  respect  to  the 
surface  normal,  the  transmission  characteristics  of  the 
dielectric  interfaces  are  monotonically  decreasing 
functions  of  incidence  angle.  Hence,  for  applications 
in  which  only  H-plane  element  pattern  or  grating  lobe 
control  is  required,  spatial  filters  are  readily 
synthesized  according  to  the  prescription  in  Figure  7, 
provided  the  slab  thicknesses  are  not  so  large  as  to 
introduce  spurious  far  out  passband  as  described  by 
Mailloux^.  However,  for  the  case  of  E-plane  control, 
or. the  more  general  case  of  providing  control  in  two 
dimensions,  the  filter  performance  is  fundamentally 
limited  by  the  Brewster  angle  associated  phenomena  as 
discussed  in  Section  2.1.  Consequently,  the  synthesis 
procedure  developed  here  has  been  specif ically  designed 
to  limit  the  spatial  extent  of  the  far  out  passband  which 
occurs  for  TM  incidence  while  simultaneously  maximizing 
transmission  in  the  passband  around  normal  incidence. 

From  the  above  considerations,  the  formal  state- 
ment of  the  synthesis  goal  is  to  minimize  the  function 
U over  a properly  defined  multidimensional  feasible 
region  in  f t } and  f t ^ } . The  function  U is  defined  as 

. o up, 

25)  U=Max{  Min  ( j T (9.,  f ) | -10  ),0) 

O<0<91 

f. <f<f 

1 n 

• 05p2 

+Max{  Max  ( |.T  (6  , f ) ; -10  ),0} 

0_<9<6 

2 s 
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where 


Maxi  a,b}  = 


b, 

a. 


for  b>a 
for  a>b 


Min  (x  (8 , f ) ) = 
0<8<81 

fl<f<fh 


minimum  value  of  x(8,f) 
in  the  two  dimensional 
space,  bounded  as  indicated 


Max  (x 
82<8<8 

fx<f<f 


(e,f ) 

s 

h 


) = 


maximum  value  of  x(8,f) 
in  the  two  dimensional 
space 


f^  and  f^  are  the  lower  and  upper  frequency  bounds, 
respectively;  8S  is  the  angle,  6S  > 02  at  which  the 
slope  of  |T(6,f)|  with  respect  to  8 changes  from  negative 
co  positive;  and  T(0,f)  is  the  voltage  transmission 

o 

coefficient  given  by  equation  (13)  for  TM  incidence. 

For  the  special  case  wherein  only  TE  polarization  is  of 
interest,  8g  = 90°  and  T(8,f)  is  taken  as  the  TE  trans- 
mission coefficient. 

The  effect  of  the  first  term  in  equation  (25) 
is  to  force  the  transmission  loss  in  the  spatial  passband 
to  be  less  than  db  over  the  frequency  band.  The 
effect  of  the  second  term  is  to  simultaneously  maximize 
the  filter  rejection  over  the  spatial  region  02<0<8S 
and  drive  9S  toward  90°. 

To  ensure  that  the  ultimate  filter  design  represents 
an  optimum,  yet  practical,  configuration,  the  feasible 
region  in  {e^}  and  {t^}  is  bounded  and  discontinuous. 

The  allowed  values  of  slab  relative  permittivity  are 
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taken  from  a table  of  available,  low  loss  microwave 

materials,  and  are  bounded  by  = 1 (vacuum)  and  a 

moderately  large  value  available  in  practice  (i  ^.  = 25 

is  a practical  upper  bound  at  X band) . The  allowed 

values  of  slab  thickness  are  continuous  in  a range  guaranteeing  a 

design  which  may  be  fabricated  using  current  manufacturing 

techniques  for  large  surfaces.  A lower  bound  of  . 020*' 

inches  is  representative,  and  an  upper  bound  of  1.0" 

is  a reasonable  engineering  limit  for  x band  designs. 

In  principle,  the  filter  synthesis  may  proceed 
from  this  point  without  further  restrictions  on  the  solu- 
tion space.  However,  as  was  discussed  in  Section  2.1, 
the  desired  filtering  characteristics  are  best  achieved 
for  permittivity  ratios  p>>l,  or  p<<l,  at  the  dielectric 
interfaces.  Such  ratios  are  obtained  by  constructing 
the  filter  with  alternating  layers  of  high  (r  1)  and 
low  (e^al)  permittivities.  Consequently,  the  synthesis 
procedure  is  designed  to  consider  only  filters  of  the 
type  shown  in  Figure  8,  which  consist  of  N dielectric 
layers  with  permittivities  {e^},  i = 1,...,N,  and  thick- 
nesses {t.},i  = 1 , . . . , N separated  by  N-l  layers  of  air 
or  low  permittivity  (tjwl)  dielectric. 

In  the  figure,  the  high  permittivity  layers  are 
shown  to  have  arbitrary  e...  It  was  found  during  the 
course  of  tnis  study  that  the  synthesis  procedure  typically 
resulted  in  a disposition  of  permittivities  which  is 
nearly  symmetric  about  the  geometric  center  of  the  filter. 

Since  this  is  not  at  all  unreasonable  in  light  of  the 
results  obtained  by  Mailloux^,  and  since  it  is 
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Figure  8 - Generic  Configurations  for  Stratfied  Dielectric 

Spatial  Filters 

desirable  to  limit  the  number  of  different  materials 

in  the  filter,  the  last  restriction  imposed  on  the 

solution  space  disposes  the  permittivities  symmetrically. 

In  broad  outline,  the  numerical  synthesis  procedure 

is  based  on  the  RAZOR  search  optimization  described  by 

(2) 

Bandler  and  MacDonald  . For  a given  incident  polarization, 
(typically  TM  for  most  applications),  frequency  and  band- 
width, desired  spatial  performance  criterion  (e.g.  Figure 
7),  and  number,  N,  of  high  permittivity  dielectric  layers, 
a modified  pattern  search  is  conducted  to  locate  local 
and/or  global  minima  of  the  function  U (equation  (25)) 
within  the  bounded  feasible  region  of  {e^},  ft^},  where 
the  are  symmetrical ly  disposed  with  j = l,...,N/2; 
and  the  t-  are  the  thickness  parameters  of  each  filter 
layer  with  i=l,...,2N-l  (thicknesses  of  both  high  and 
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low  permittivity  layers  are  considered) . An  estimate 
{t°}  is  used  to  initiate  the  search.  If  the 
pattern  search  does  not  locate  a minimum  of  U which 
satisfies  a predetermined  measure  of  success,  f , that  is 


26) 


U > 


f 


c 


a new  pattern  search  is  initiated  by  randomly  perturbing 
the  estimate  {e°} , {t°}  with  large  steps  along  the 
multidimensional  coordinate  axes.  After  a set  number 
of  evaluations  of  U,  or  after  a set  number  of  failures, 
as  defined  by  equation  (26),  the  procedure  concludes, 
identifies  the  cause  of  termination,  and  provides  the 
best  estimate  { e ^ } , { } as  defined  by  the  minimum  value 
of  U-f  . If  a success  is  recorded,  i.e.,  U<f  , the 
procedure  terminates  and  provides  the  solution  set  { e . } , 
3=1 , ... , N/2 , and  {t^},  i=l , . . . , 2N-1 . 

A detailed  account  of  the  procedure  is  given  in 
Appendix  A. 
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3.  NUMERICAL  RESULTS 


In  this  section,  numerical  results  are  presented. 

In  the  first  part,  results  are  given  which  illustrate 
the  manner  in  which  the  synthesis  technique  approaches 
the  optimum  solution.  In  the  second  part,  a spatial 
filter  design  is  given  which  provides  far  out  sidelobe 
and  quantization  lobe  control  for  a large  mechanically 
scanned  planar  array. 

3. 1 Synthesis  of  Seven  Layer  Filters 

Figures  9 through  16  show  typical  results  obtained 
for  seven  layer  filter  at  various  stages  of  the  synthesis 
procedure  developed  in  this  study.  The  synthesis  is 
carried  out  only  for  TM  polarization  since,  as  demonstrated 
above,  filteiing  of  incident  TE  waves  does  not  present 
a significant  design  problem.  In  obtaining  these  results, 
slab  relative  permittivities  were  allowed  to  take  on 
any  value  in  the  range  l<e.^<25,  and  thus  the  configurations 
are  not  practical  in  the  sense  defined  previously.  The 
synthesized  designs  are  given  only  to  illustrate  the 
result  of  repeated  application  of  the  procedure  to  obtain 
the  desired  filtering  characteristic.  The  performance 
goal  is  indicated  on  each  figure.  In  the  spatial  range, 

0- 6<10°,  the  transmission  loss  is  less  than  ldb.  In 
the  range,  25°<6,  the  desired  rejection  is  10  db  or 
greater.  The  frequency  band  is  4%  about  11.8  GHz. 

The  initial  estimates  {i.j),j  = l,2,  and  { t ^ , 

1- l,..., 7 were  taken  from  the  filter  design  given  by 
Mailloux^.  The  outer  high  permittivity  slabs  have 


24 


5INB 


LEGENI 

Q FREQ  = .“E  q,  FREQ  =r 

a FREQ  r:  1.Q2 

Figure  9 - TM  Power  Transmission  of  Seven  Layer  Spatial 
Filter  (from  Mailloux^) 
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□ FREQ  = .5»B  FREQ  = X.DD 

a FREQ  = X.D2  . 


Figure  12  - TE  Power  Transmission  of  Seven  Layer  Spatial  Filter 
e .=3. 95, £,=14. 91, (t. }={ .142, .440, .046, . 500, .057, .965, .129} 
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Figure  13  - TM  Power  Transmission  of  Seven  Layer  Spatial  Filte 

6^3.93,^  =14.93,  (t  } = {.  139,  .440,  . 046,  .500,  . 057,  . 965,  . 129  • 
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Figure  14  - TE  Power  Transmission  of  Seven  Layer  Spatial  Filter 
e =3.9 3, e =14. 93, it  }=  l . 139 , . 4 40 , . 046 , . 500 , . 057 , . 965 , . 129 
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Figure  15  - TM  Power  Transmission  of  Seven  Layer  Spatial  Filter 
e ,= 7. 27, 6-  = 18. 75, {t.}={. 114, .449, .051, .483, .050, .947, .074} 
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q FREQ  = .9S  * FREQ  = X.DD 

<♦,  FREQ  = 1.Q2 

Figure  16  - TE  Power  Transmission  of  Seven  Layer  Spatial  Filter 
e , = 7. 27 , e,=18. 75, (t . } = {.  114, . 44  9, . 051, . 4 83, . 050,  . 94  7,  . 074} 


e^=3.08,  and  the  inner  high  permittivity  slabs  have 
e2=15.14.  The  layer  thicknesses  are  .25//ei  inches  in 
the  high  permittivity  slabs,  and  1.0"  in  the  intervening 
regions  which  are  assumed  to  be  free  space.  The  per- 
formance of  ti»is  filter  for  TM  and  TE  incident  polariza- 
tions is  shown  in  Figures  9 and  10,  respectively.  For 
TM  and  TE  incidence,  tne  filter  provides  the  required 
low  transmission  loss  in  the  spatial  passband  and  over 
the  entire  frequency  band.  However,  the  spatial  stop- 
band  is  very  narrow  for  TM  incidence  and  a spurious 
farout  passband  appears  for  both  polarizations  near  60° 
incidence  angle.  The  goal  of  the  numerical  optimization 
technique  is  to  eliminate  the  spurious  passbands  for  both 
polarizations  and  extend  the  spatial  stopband  for  TM 
incidence . 

Figures  11  and  12  show  the  filter  response  after 
200  evaluations  of  the  function  U.  It  is  evident 
unat  the  synthesis  procedure  is  obtaining  greater  stop- 
band  extent  at  the  expense  of  transmission  loss  in  the 
passband.  For  both  polarizations,  the  spurious  passband 
has  been  significantly  reduced.  The  thickness  distribution 
and  relative  permittivity  have  been  sign  icantly  altered. 

After  400  evaluations  of  the  function  U,  the 
synthesized  filter  response  is  as  shown  in  Figures  13 
and  14.  Little  change  has  occurred  with  respect  to  the 
previous  result,  and  it  is  clear  that  a continued  search 
in  the  immediate  vicinity  will  not  obtain  the  desired 
result.  Using  the  current  best  estimate  as  the  starting 
point,  a large  step  size  is  introduced  for  the  pattern 
search  to  drive  the  result  away  from  the  local  minimum. 

And,  after  an  additional  200  function  evaluations,  the 
performance  shown  in  Figures  15  and  16  is  obtained, 


33 


which  satisfies  the  design  goals.  In  particular,  the 
extent  of  the  stopband  has  been  increased  from  the  original 
16°,  beginning  at  6-32°,  to  31°  beginning  at  8=26°. 

Although  significant  transmission  loss  has  been  intro- 
duced in  the  passband,  it  may  be  removed  by  further 
optimization  in  the  vicinity  of  the  current  best  esti- 
mate. It  is  particularly  interesting  to  note  that  the 
farout  TM  incidence  passband  peak,  has  been  moved  well 
outside  the  first  layer  Brewster  angle  and  the  width 
at  the  5 db  points  has  been  reduced  from  30°  for  the 
original  design  to  15°. 

3 . 2 An  Eleven  Layer  Spatial  Filter  for  Use  with  a Large 
Mechanically  Scanned  Planar  Array 

In  this  section  a practical  filter  design  for  use 
with  a large  planar  non-scanning  array  is  presented. 

Since  the  permittivity  of  most  microwave  materials  is 
constant  from  S band  to  X band,  the  thickness ' parameters 
determined  by  the  synthesis  procedure  are  given  normalized 
to  the  free  space  wavelength  at  the  center  frequency  of 
the  operating  band,  and  the  design  may  be  considered 
universal  for  this  frequency  range. 

The  design  goal  is  to  provide  a minimum  farout 
(8>25°)  sidelobe  and  quantization  looe  reduction  of  10 
db  for  a large,  mechanically  scanned,  two-dimensional 
planar  array  of  vertically  polarized  rectangular  apertures 
over  a 4%  frequency  band.  The  array  possesses  full 
search  track  capability  and  requires  a 10°  spatial 
passband  with  maximum  transmission  loss  of  . ldb.  The 
required  coverage  sector  is  360°  azimuthally,  and  -5° 
to  50°  vertically.  An  SOJ  threat  is  postulated  to  be 
uniformly  distributed  in  azimuth  within  20°  of  the 
horizon. 


34 


To  provide  the  required  rejection  over  the  full 
■coverage  sector,  the  filter  design  must  limit  the  Brewster 
angle  associated  passband  to  angles  greater  than  55°. 
Practical  filters  providing  a 30°  E-plane  stopband  may 
be  synthesized  using  eleven  or  more  layers. 

Figure  17  shows  an  illustration  of  an  eleven 
layer  filter  Wi.ich  satisfies  the  design  goals.  The 
filter  is  constructed  using  four  well  known  dielectric 
materials:  stycast  Hi-K  (£^  = £^=5)  ' an  t:merso11 

Cumming  loaded  cross-linked  polystrene;  Trans-Tech  DA-9 
Alumina  (e.  q=9 . 5 ) ; stycast  Hi-K  500F  loaded  uhermoset 

hydrocarLon;  and  (e _=e „=£=£=£ . n=l . 02 ) . The  layer 
thicknesses  ‘are  given  normalized  to  free  space  wavelength 
at  center  frequency. 

The  filter  performance  is  shown  in  Figures  18  and 
19  for  E and  H plane  incidence,  respectively.  In  the 
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Figure  17  - Cross  Section  of  Filter  Synthesized  by 
Numerical  Optimization 
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Figure  18  - TM  Power  Transmission  of  Eleven  Layer  Spatial  Filter 
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Figure  19  - TE  Power  Transmission  of  Eleven  Layer  Spatial  Filter 


spatial  passband,  the  transmission  loss  is  less  than  .15 
db  over  the  4%  frequency  band.  The  E-plane  spatial  stop- 
band  extends  from  approximately  18°  to  61°  at  low  freq- 
uency and  from  26°  to  62°  at  the  high  frequency  and 
rejection  is  typically  14  db  or  greater  from  30°  to 
55°.  Excellent  stopband  rejection  is  obtained  in  the 
H-plane. 

Figure  20  shows  a center  frequency  contour  plot  in 
sine  space  of  filter  transmission  in  db  for  plane  waves 
incident  at  angle  6,0  with  respect  to  the  coordinate 
system  of  Figure  2-1.  The  plane  wave  generator  is 
assumed  to  be  an  X-directed  magnetic  current  on  the 
array.  In  the  passband,  transmission  is  seen  to  be 
roughly  independent  of  6 and  0,  and  consequently  the 
filter  will  not  distort  the  main  beam.  In  the  stopband, 


Figure  20  - Transmission  of  Fields  Generated  by  X-Directed 
Magnetic  Currents  on  the  Array 
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the  rejection  is  greater  than  10  db  except  in  the  vicinity  or 
the  near  in  stopband  edge  for  0=90°.  Over  most  of  the  step- 
band,  the  rejection  considerably  exceeds  14  db . The 
Brewster  angle  associated  passband  does  not  extend  to 
angles,  0,  below  33°,  thereoy  reducing  the  SOJ  threat/ 

The  isolation  between  nominal  (or  transmitted) 
field  polarization,  and  the  cross  polarization  generated 
by  the  filter  is  shown  in  Figure  21.  As  expected,  the 
peak  cross  polarized  signal  is  generated  in  the  vicinity 
of  0=45°,  0=20°,  and  in  the  farout  passband  region. 


RKT  A VA»A?.!F  fftPY 

-.L  Lui  I 


• 2 • 2n  i 

sm  a * sir  B - 1 


On 


Figure  21  - Isolation  of  Cross  Polarized  Signal  Generated  by 

Layered  Medium 


4.  CONCLUSIONS 


The  stratified  dielectric  slab  spatial  filter  syn- 
thesis technique  given  here  results  in  practical  filter 
designs  which  provide  considerable  element  pattern  control 
for  large  mechanically  scanning  and  limited  scanning 
arrays.  For  a hypothetical  mechanically  scanned  two 
dimensional  antenna,  a filter  design  has  been  presented 
which  reduces  farout  sidelobes  and  quantization  lobes 
by  more  than  10  db  in  a region  extending  from  26°  off 
the  antenna  normal  to  61°. 
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APPENDIX  A 

A MODIFIED  RAZOR  OPTIMIZATION  PROCEDURE 
FOR  SPATIAL  FILTER  SYNTHESIS 

I 

i 

A . 1 Program  Abstract 

The  optimization  program  is  a federation  of  three 
modules  named  FILTER,  RXS,  and  RZS  which  will  be  described 
as  functional  entities.  FILTER  is  the  main  program 
segment  containing  in  addition,  a subroutine  named  U 
which  computes  the  objective  function  to  be  optimized, 
and  DISPLA,  a utility  subroutine  which  surpresses  the 
automatic  carriage  return  and  line  feed  after  printing 
a prompting  message. 

The  main  program  permits  the  user  to  input  all 
pertinent  parameters  such  as:  the  number  of  slabs  of 

the  spatial  filter,  the  initial  dielectric  constants 
and  thicknesses,  tolerance  ranges  over  which  the  opti- 
mization is  to  seek  its  improved  values,  the  frequency 
ranges  over  which  the  filter  is  to  operate,  the  sine 
theta  and  transmission  coefficient  values  which  character- 
ize the  pass/stop  bands  and  the  performance  in  db,  as 
well  as  the  parameters  ALPHA,  LIMIT,  EPSMIN,  RHO,  ETA, 
KAPPA,  DELTA,  FEPS  which  govern  the  conduct  of  the  opti- 
mization process  itself. 

Sufficient  flexibility  has  been  incorporated  into 
the  main  program  to  allow  the  user  to  select  a TRACE 
option  which  permits  witnessing  the  deliberations  of  the 
optimization  by  a printout  of  all  perturbated  components 
and  their  effect  on  the  objective  function.  In  addition, 
the  program  has  been  structured  so  that  only  a minor 
alteration  by  the  user,  if  he  so  desires,  will  enable 
him  to  treat  special  filter  configurations  like  symmetric 
filters,  quarter  wave  length  layer  filters  and  so  on. 
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The  second  module  RXS  consists  of  a single  sub- 
routine RXTX  which  returns  the  complex-valued  transmission 
coefficient  RT  used  only  by  the  function  U.  RT  results 
from  the  matrix  multiplication  of  junction  and  line  length 
wave  transmission  matrices,  and  as  this  coefficient  is 
computed  many  times  by  the  function  U (even  for  a single 
perturbation) , great  care  has  been  exercised  in  the 
construction  of  the  algorithm  for  the  total  wave  trans- 
mission matrix  computation. 

The  last  module  RZS  consists  of  six  subroutines 
named  RAZOR,  FINISH,  PATSER,  BOUND,  EXPLR , PATMV , with 
FINISH  and  BOUND  serving  in  ancillary  capacities  of 
determining  whether  the  perturbations  have  produced  a 
minimization  of  the  objective  function  and  secondly 
ensuring  that  no  perturbation  is  allowed  to  exceed  its 
prescribed  range.  RAZOR  is  the  subroutine  which  in  con- 
junction with  PATSER,  PATMV,  EXPLR  executes  the  optimiza- 

(2) 

tion  technique  described  by  of  Bandler  and  Macdonald 
In  its  broad  outline  the  technique  is  to  conduct  a pattern 
search  (PATSER)  around  a base  point  by  perturbing  in  turn 

O 

each  of  the  components  of  the  base  point.  If  the  search 
exposes  a new  point  for  which  the  objective  function  has 
been  minimized  to  within  a tolerance  specified  by  FEPS, 
RAZOR  returns  this  new  point  to  the  main  program  and 
terminates  its  activities.  Otherwise,  a point  in  the 
vicinity  is  randomly  selected  as  a provisional  base 
point  around  which  explorations  are  undertaken  as  before. 
If  there  is  no  improvement,  the  process  is  repeated  for 
as  many  times  as  specified  by  KAPPA  with  a closer 
randomly  selected  provisional  base  point. 
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A. 2 Program  Description 


The  main  program  TEST  1 calls  for  user  input  cf 
NS,  the  number  of  non-air  dielectrics,  the  associated 
dielectric  constants,  slab  thicknesses,  and  the  high  and 
low  constraints  on  the  perturbation  vector  PHO.  In 
addition,  the  frequency  range  over  which  the  filter  is 
to  perform  is  input  as  are  the  sine  theta  and  transmission 
coefficient  values  characterizing  the  pass/stop  bands  and 
performance  in  db . 

Thus,  for  a spatial  filter  with  six  non-air 
dielectrics  configured  symmetrically,  6 would  be  entered 
for  NS  in  response  to  the  prompt  NR  SLABS=?.  Also 
3 values  ( IP= (NS+1 ) /2 ) of  non-air  dielectric  constants 
and  11  values  (NSLAB1=2*NS-1 ) of  air  and  non-air  thicknesses 
would  be  input.  The  dimensionality  of  the  perturbation 
vector  PHO  would  then  be  14  (K=IP+NSLAB1 ) , and  whose 
initial  contents  would  be  the  3 dielectric  constants  and 
11  thickuesses.  The  next  two  input  parameters  would  be 
two  sets  “of  K (here,  14)  values  for  low  and  high  pertur- 
bation constraints  to  be  placed  on  the  permissible  ranges 
of  the  PHO  vectors  components.  Typically,  a low  of  1 
and  a high  of  25  have  been  placed  on  the  PHO  components 
associated  with  dielctric  constants,  a low  of  .02  and 
high  of  1.0  for  those  components  associated  with  non-air 
dielectric  thicknesses,  and  a low  of  .1  and  a high  of  2.0 
for  air  thicknesses.!  In  response  to  the  prompt  for  freq- 
uency range,  the  user  enters  3 values  for  low  frequency, 
high  frequency,  frequency  increment  normalized  to  center 
frequency;  for  example  .48,  1.02,  .02  would  cover  three 
frequencies.  The  next  six  values  to  be  entered  are 
concerned  with  the  filter's  performance.  The  value 
specified  for  DB1  determines  the  low  cutoff  figure  for 
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transmission  coefficient  in  the  passband  and  the  values 
entered  for  SN1  and  SNINC1  in  response  to  the  prompt 
STH1  and  I ICR=?  define  the  spatial  limit  of  the  passband 
and  the  sequence  of  points  within  the  passband  which  are 
to  be  evaluated.  Similar  remarks  apply  to  DD2,  SN2, 

SNINC2  whic.i  specify  the  stopband.  Values  of  DBl=-.l, 

SNl= .175,  SNINC1=. 0175 , DB2=-10.0,  SN2=.42,  SNINC2=.02 
are  typica..  The  program  next  inputs  ALPHA,  LIMIT,  EPSMIN, 
RHO,  ETA,  DAPPA,  DELTA,  and  FEPS.  FEPS  is  the  tolerance 
within  which  our  objective  function  is  to  be  minimized 
and  is  typically  selected  as  .01.  RAZOR  conducts  a pattern 
search  (PATSER)  around  the  initial  base  point  by  perturb- 
ing in  turn  each  of  the  components  of  the  base  point.  The 
actual  absolute  step  sizes  used  in  PATSER  are  SDEL=DELTA* 

R ( I ) *S ( I ) /100  where  R(I)  = 100* (HIGH ( I ) -LOW ( I ) ) is  the 
parameter  range  of  permissible  values  in  the  Ith  component, 
and  S ( I ) has  values  +1  or  -1  to  establish  directionality 
of  the  perturbation. 

Since  dielectric  constants  are  considered  in  the 
range  l-*-25  and  thicknesses  .02-+ 1 for  non-air  dielectrics 
and  . 1-+2  for  air  dielectrics,  we  observe  that  selecting 
an  initial  DELTA= . 1 yields  step  sizes  of  2.4  for  dielectric 
constants  with  .008  for  non-air  thicknesses  and  .19  for 
air  thicknesses. 

EPSLON  and  DELTA  are  used  as  relative  measures  of 
step  size  reductions,  with  DELTA<EPSLON  signifying 
? that  a pattern  search,  PATSER,  would  be  useless  at  this 

point.  Initially  EPSLON=EPSMIN* (ETA**KAPPA) 

where  EPSMIN  = .001  (a  starting  value) 

ETA  =2  (a  halving  factor) 

KAPPA  = 3 (max  of  3 random  moves) . 
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EPSMIN  = .001  was  selected  so  that  EPSLON  = .001  *2i  - 
thereby  rendering  the  initial  DELTA  = . 1 larger  than  EPSLON 
and  thus  guaranteeing  a PATSER  the  first  time  through. 

ALPHA  is  a parameter  within  PATSER  which  reduces 
the  size  of  DELTA  in  the  event  initial  explorations  give 
no  improvement  in  reducing  the  value  of  the  objective 
function.  A va.ue  of  .25  gives  significant  reduction  of 
DELTA,  but  does  not  reduce  coo  rapidly  to  obtain  a 
solution . 

When  all  explorations  about  the  base  point  have 
been  unsuccessful,  RAZOR  perturbs  each  component  of  the 
base  point  by  RHO*RANDOM*EPSLON/1QO  in  order  to  obtain 
a new  provisional  base  point  about  which  the  process 
can  be  repeated.  RANDOM  provides  random  values  between 
-1  and  +1  with  RHO  selected  as  400.  EPSLON  is  then  halved 
by  ETA  and  new  explorations  are  begun. 

Since  an  enormous  number  of  functional  evaluations 
must  be  performed  not  only  in  exploratory  moves  but  in 
evaluatory  criteria,  LIMIT  (the  maximum  number  of  function- 
al values  to  be  computed  in  the  entire  optimization 
program)  has  been  set  to  a large  value,  typically  on 
the  order  of  1000.  The  objective  function  U is 
summarized  in  Figure  7 and  in  equation  25.  Its 
actual  implementation  is  performed  in  two  steps.  For 
each  fixed  sine  theta  (sine  theta  = 0 through  SN1  in 
steps  of  SNINC1 ) cycle  through  the  range  of  frequencies 
FLO  through  FHI  in  steps  of  FINCR  constantly  determining 
UMN=CABS (RT)  - OFFSET  where  OFFSET  = 10.0  + (.05*DBl) 
and  RT  is  the  complex  valued  transmission  coefficient 
returned  from  a call  to  subroutine  RXTX.  The  negative 
of  the  minimum  value  of  these  UMN  values  is  kept  as  the 
first  part  of  the  objective  function. 
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The  second  part  is  determined  in  the  following  manner. 
Beginning  with  a fixed  sine  theta  value  of  SN2,  cycle 
through  the  range  oi  1 1 equenc  i es  FLO  through  Kill  in  stops 
of  FINCH  constant  ly  dc  1 oniii  n i ng  UMX  ('All.'!  (RT)  - OFFSET  whoto 
OFFSET  = 10.0  t (.05*DB2).  Continue  this  process  with 
the  next  fixed  sine  theta  value  incremented  by  SN1NC2 
over  its  previous  value  only  for  as  long  as  th^  center 
frequencies  CABS (RT)  values  are  declining. 

A. 3 Program  Flow  of  Principle  Programs  and  Subroutines 

To  give  a more  complete  description  of  the  optimiza- 
tion procedure,  the  following  paragraphs  simultaneously 
outline  program  flow  and  comment  on  the  implementation  of 
the  principle  programs  and  subroutines.  RAZOR (0°,  emin»P» 

n, kappa)  is  the  principal  procedure  in  the  optimization 

„o 

process.  0 , cmin'P'n»  kappa  are  FORMAL  parameters 
through  which  the  initial  values  are  acquired  from  the 
calling  program.  Within  the  algorithm  0,  U0  are  local 
variables,  k is  the  dimensionality  of  the  0 vector,  and 
6 is  common  (i.e.  global)  to  RAZOR  and  has  been  initialized 
by  the  main  program  to  a value  6>e  *r)^caPPa  to  ensure 
bypassing  the  IF  6<e  criterion  in  PATSER  whica  RAZOR  invokes 
at  the  outset  of  its  deliberations. 

Comment:  initialize  the  local  quantity  U0°  to  reflect 

the  starting  value.  Set  up  the  initial  direc- 
tions for  the  explorations.  Initialize  the 
starting  t and  then  invoke  a pattern  search. 

If  a point  is  uncovered  with  a functional 
value  within  the  FEPS  tolerance  specified, 
then  terminate  the  procedure  otherwise  renew 
the  search  (kappa  number  of  times  if  necessary) 
with  an  arbitrary  point  randomly  selected  in 
the  vicinity  of  0°  and  with  a reduced  e value; 
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U0°:  = U(0°);  for  i:=l  step  1 until  < DO  S^:=l; 

e:=emm*nkdPPa;  PATSER  {P°'  u^°); 

If  finish  (Uj0° ) then  go  to  FIN; 

for  j:  = l step  1 until  kappa  DO 

Begin  comment  initialize  a local  variable  U0°  to  the  best 
functional  value  and  a local  vector  0 to 
contain  the  coordinates  of  a randomly 
selected  point  in  the  vicinity  of  0°.  Another 
local  variable  U0  is  initialized  to  contain 
the  functional  value  at  the  random  point, 
e and  6 are  reduced  and  a pattern  search 
is  initiated; 


U0°:=U(0°) ; 

for  i:=l  step  1 until  k DO  0^:=0°  +P  *RANDOM*e; 

c:=e/n;  6 : = | 0-0° ! /k ; 

U0 : =U (0 ) ; PATSER  (0,U0); 

If  U0<U0° 

Then  LI:  begin  comment  set  U0°  to  the  returned  improved 

value  U0.  Establish  0 as  the 
direction  of  improvement  and  0° 
to  the  returned  improved  location 
f'-o 

U0° : =U0 ; 0:=0-0°;  0°:=0 
End 
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Else 


Begin  Comment  reset.  1)0  to  the  starting  functional  value 
U0°  over  which  there  was  no  improvement 
and  set  8 in  the  opposite  random  direction; 

U0:=UJ3°;  8:=0°-0 

End; 

Comment  set  local  0 as  the  new  extrapolated  point  from 
0°  in  the  8 direction  and  6 equal  to  the  square 
of  the  magnitude  of  8.  Invoke  a pattern  move; 

0 =0°+8 ; 6 : = 1 8 | 2 ; PATMV  (U0,0,0°); 

Comment  if  the  returned  functional  value  is  an  improve- 
ment then  return  to  LI  for  further  refinement 
until  none  better  is  found.  Then  check  if  the 
returned  refinement  is  within  the  FEPS  tolerance 
which  causes  a termination  of  the  RAZOR  optimiza- 
tion. Failing  this,  the  next  randofo  point  is 
selected  even  closer  to  0° ; 

If  U0<U0°  then  go  to  LI; 

If  finish  (U0° ) then  go  to  fin 

End; 

Write  ("no  convergence"); 

Fin:  End 


PATSER  (0° , U0° ) is  the  chief  searching  procedure.  0° 
and  U0°  are  formal  parameters  and  0,U0  are  local  variables 
in  the  following  algorithm. 


A-8 


IF  6<e  then  go  to  fin; 

Comment  initialize  the  local  variables  0,U0  to  the 
corresponding  values  transferred  through 
FORMAL  parameters  0°,  U05.  Then  invoke 
EXPLR  to  perform  initial  explorations. 

0,U0  will  be  changed  hopefully  to  an 
improved  point  and  its  functional  value; 

U0:=UjZf°;  0:=0°;  EXPLR  (U0,0) ; 

If  U0>U0° 

Then  begin  comment  there  is  no  improvement  in  any 

of  the  explorations.  Reduce 
size  of  6 and  try  again; 

6:=a*6;  go  to  LI 

End; 

Comment  an  improvement  was  detected.  Change  the 
formal  parameters  0°,U0°  to  reflect  the 
improved  location  and  its  functional 
value.  Establish  into  0 the  general 
direction  of  improvement.  Set  the  local 
quantity  0 to  contain  the  new  extrapolated 
value  and  6 to  contain  the  square  of  the 
magnitude  of  0.  Invoke  a pattern  move 
(PATMV)  and  return  to  L2  if  there  is  an 
improvement,  otherwise  begin  anew  by 
returning  to  LI; 


L2 : U0°:=U0;  6: =0-0°;  0°:=0;  0:=0°t0; 

5 : = 1 0 i 2 ; 

PATMV  ( U0 , 0 , 0° ) ; 

If  U0- U0°  then  go  Lo  L2  else  go  to  LI; 

Fin:  End 

PATMV  (U0° ,0,0^)  where  0°  is  the  vector  describing  the 
best  location  so  far  and  0 is  the  extrapolated  location 
around  which  explorations  are  to  be  performed  within 
PATMV,  is  best  described  graphically  in  terms  of  a two 
dimensional  0.  If  the  explorations  depicted  below  un- 
cover a better  point  than  0,  then  that  point  is  returned 
in  place  of  0 and  the  functional  value  at  the  uncovered 
point  is  returned  to  U0.  Within  the  algorithm  K is  the 
dimensionality  of  the  0 vector,  while  6 and  the  vector 
0 are  common  (i.e.  global)  to  PATMV.  The  functional 
value  of  the  original  extrapolated  point  0 is  computed 
immediately  upon  entry  into  PATMV  and  is  retained  in  a 
local  variable  named  UPHI.  It  is  against  this  local 
quantity  that  functional  values  of  exploratory  test 
points  (about  0)  returned  by  EXPLR  are  compared  to 
determine  if  improvement  has  occurred. 

PATMV  first  performs  explorations  (using  EXPLR) 
with  a new  5=/6/<  and  returns  into  0 any  of  the  possible 
locations  indicated  and  into  U0  the  corresponding 
functional  value.  Note  that  0 is  the  direction  estab- 
lished earlier  in  the  process  and  reflects  the  general 
direction  of  improvement  from  some  former  0°  location 
to  the  present  0°  location  specified  upon  entry  into 
PATMV. 
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Figure  A-l  - Initial  Exploration  in  Subprogram  PATMV 


If  the  exploratory  moves  about  j2  fail  to  uncover  a better 

point  or  in  the  case  there  is  an  improvement  but  there 

-S  no  significant  change  in  each  component  i.e.  for 

every  i |0.-0o)<lO  6*j0°  , 

1 i 1 1 1 l 1 , then  the  step  size  is  halved 

to  1/2 /$/<'  and  explorations  are  conducted  about  a new 

extrapolated  location  half  as  far  away  from  0°  in  the 

0 direction  as  shown. 


? ? 


Figure  A-2  - Explorations  About  New  Extrapolated  Vector 

0 

• $ = j + $ 
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Finally,  if  there  is  still  no  improvement  or  in  the  case 
there  is  an  improvement  but  there  is  no  significant 
change  in  every  component  J2T  in  the  aforementioned  sense, 
then  the  halved  step  size  1/2/6/k  is  used  again  but 
explorations  are  conducted  about  a new  extrapolated 
location  half  as  far  away  from  0°  in  the  negative  6 direction 
as  shown. 


Figure  A-3  - Explorations  About  New  Extrapolated  Base  Point  in 

-6  Direction 


EXPLR  (U0°,0) 

U0°  and  0 are  formal  parameters.  Within  EXPLR  U0  is 
a local  quantity.  EXPLR  will  change  U0°  and  0 to  an 
improved  functional  value  and  an  improved  point,  re- 
spectively. < is  the  dimensionality  of  the  0 vector. 
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In  graphic  terms  of  a two  dimensional  0 , EXPLR  will 
return  into  0 any  of  the  possible  positions  indicated, 
and  into  U0°  the  corresponding  functional  value. 
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Figure  A-4  - Two  Dimensional  Graphic  Representation  of  EXPLR 

The  center  point  is  the  location  of  0 upon  entry  into 
I EXPLR.  In  the  event  that  EXPLR  is  unable  to  find  an 

improvement  in  U0°  at  any  of  the  locations  sdel  away 
from  the  center  point,  upon  exit  from  EXPLR,  the  original 
functional  value  U0°  and  the  original  0 location  remain 
unchanged. 

FINISH  (U0°)  is  the  logical  procedure  which  furnishes  a 
value  of  true  if  the  functional  value  contained  in  its 
formal  parameter  is  within  the  tolerance  FEPS  initialized 
by  the  mam  calling  program. 
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APPENDIX  B 

MICROWAVE  DIELECTRICS  FOR  APPLICATION  IN 
STRATIFIED  DIELECTRIC  SLAB  SPATIAL  FILTERS 

o 

The  table  given  in  this  appendix  summarizes  the 
availability  of  low  loss  microwave  dielectrics  suitable 
for  use  in  stratified  dielectric  slab  spatial  filters. 
The  data  presented  was  obtained  over  the  period  of 
1 Apr  1976  to  1 Oct  1976,  and  as  such  represents  an  up 
to  date  compilation  of  available  data.  Where  possible, 
information  relative  to  physical  properties  is  included. 

All  materials  in  the  table  have  constant  or  nearly 
constant  electrical  properties  from  S band  to  X band. 

The  loss  tangent  of  all  materials  is  less  than  .002. 
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APPENDIX  C 
PROGRAM  LISTINGS 


In  this  appendix,  program  listings  are  given  for 
all  computer  programs  written  in  the  course  of  this 
study.  The  listings  are  self-explanatory.  Information 
relative  to  the  optimization  programs  and  subprograms 
is  given  in  Appendix  A. 
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PROC-Ra^  TP  ST  1 (INPUT,  TAPE  j » OUTPUT  ♦ T APE6»0UT  PUT ) 

C»,,6,T>-I3  is  Tm p VERY  FIRST  WORKING  OPTIMIZATION  PROGRAM 
C # r . • » I T IS  USED  IN  CONJUNCTION  WjTh  RxS  ♦ RZS3 

COMMON  /RZR/  K. DELTA.  ALPHA, LIMIT, N»N2tN3,EPSLON, THFTAM5)  , 

• St  15} .FEPSt TRAcFt LOW ( 15) .HIGH (15) .RANGE (15) 
LOGICAL  TRACE, ODD 

RFAl.  lO* 

COMMON  /UFCT/  D(13I, FPSCUI, NS. N9(  abs.tr. odd. flo. fhi,ftncr, 

* 081 * SNI .SNINC 1 .DP2.SN2.SNINC2 
DIMENSION  M3G(3) .mSg2(?) .mSg3(?) * pHO (1 5 ) * D I ElC T ( i 0 ) 

haTA  V-SG/10HE3SMIN.RHC  . IOh.ETa. KAPPA  , 7H, DELTA*  / 

DATA  vSg?/10hF»E3  LO**I  , Bh.INcRb/ 
c..,. .INPUT  NS.  the  number  OF  NON. AIR  DIELECTRICS 

C NSLARS  IS  ThF  TOTAL  NUMBER  OF  Al»  X NON. AIR  DIELECTRICS 

C.....IP  IS  The  NUMBER  OF  DIFFERENT  NON. AIR  DIELECTRICS 

C.....K  IS  The  DIMENSION  OF  PHO  VECTOR  .TOTAL  NR  OF  PERTIIR6A  Y I OvJ,* 

1 CALL  DISPLAC9HNR  SLABS*. 1.2) 

RE'AD(l.F)  NS 

T F ( N S . LE . 01  STOP 

NSIA8S«2*N8 

NSLARl.NSL A8S-1 

IP»CNSM)/2 

K*IP*NSLAB1 

C... ..INPUT  DIELCTCI)  1 = 1 TO  IP  X INITIALIZE  EPS  VECTOR 
ENCODc(20.26R.MSG3)  IP 

289  FORMAT ( * EPS(1  TO  *,I2.*)«*) 

Call  DISPLA(msG3.2.i) 

R£  A D ( 1 * * 5 (DIELCT(T) .1*1 »IP) 
odd., TRUE. 

I F ( MOD ( NS . 2 ) .EG.  01  ODD*. False, 
do  78  T*l  , IP 
I X*2*  I 

EPS(Tx-l 5*1.0 
FPS ( I X ) *D I ElC  T ( 1 5 

IFfODO  .AND.  ( T , F G , I P ) 5 GO  TO  7* 

TX«2*CIP*I  5 
FPS(ix-n*i  .o 
lv*NS/2-M-I 

FPS(IX)aDTELCT(IY; 

7b  CONTINUE 

C... ..INPuT  DOT)  1*1  TO  N s l.  A R 1 

F-NC0DE(20.290.MSG3)  NSL  AR  1 

290  FORMAT!*  0(1  TO  *, 12. *5**5 
CALL  DISPLA (MSG3.2. l 5 
RFADd.*)  (0(1). T*1  ,N8L  ARi  ) 

C INPuT  CONSTRAINTS  l O^.MIOH.RANGf;  ENTIRE  PHO  VECTOR 

C INPUT  LOw(T)  1*1  TO  K 

ENCOOE(20.?91 .MSG35  K 

291  FORMAT ( * lOW(i  TO  *,T2.*5**5 
CALL  0TSPLA(MSG3.2.l 1 

RE  AO  ( 1 . * ) (LOw ( I ) . T s 1 . k ) 

C INPUT  h I G h ( 1 5 1*1  To  k 

encode (20*292. msg3)  k 

292  FORMAT (*hIGh( i TO  *,;2.*)**) 

CALL  0! SPL A (MSG3.2 • 1 5 

RE  AO ( 1 . * ) ( M T GH ( I ) , i * 1 « k ) 

DO  R«  T * 1 . K 

8«  RANGE  (T)*100.0*(hI6h(I)»LOMI(T)) 

C INITIAL  I7E  PhD  VECTo»  TO  INptlT  DTELECTRICS  X THICKNESSES 

DO  79  1=1. IP 
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79 


80 

C, 


c 


c 


c 


c 


V 


c 


c 


c « « 
81 

390 


8? 

C. 

83 


39 1 


C 

C 


PHO(I)*DIELCT(I3 
00  80  T*1»N3U91 
T X»  I P4>  I 
PHO(IX)SDCI) 

... INPUT  PREOLO^iFREOHi .FREOINCP 
CALL  01 SPLA  CmSG2 • 2 • J 3 
RE  AD ( 1 » * 3 FLOi PhI ,PiNCP 

...INPUT  OBI.. 

CALL  OISPLA  («M08U»  l . 1) 

RE  *D ( 1 » * 3 OBI 

...INPUT  SIN(THETAl)  X 8IN(THETAlMNrREMENT 
CALL  OISPLAdOHSTHlxINCRPf  1«1) 

REAOClt*)  SM.SNINCl 

...INPUT  OB?,.,,,,, 

CALL  DI$BLAUHDBa«*i*n 
RE  AO ( i * * 3 OB? 

...INPUT  SIN(ThETA?3  X S I N ( ThFt A? ) I Nr RE* ENT 
CALL  DISPLA(10HSTH2XINCR**1»1) 

RE  AD ( 1 * * 3 SN?»8NINC? 

...INPUT  ALPHA  X MAX  LImIT  OF  U FUNCTION  EVALUATIONS 
CALL  0ISPlA(9HALFA»lMT»i!*1  3 
RE  AD ( 1 » * 3 ALPHA*LIMlT 
, . . INPUT  EPSMIN. RhO#ETa»KAPPA »dFLTa 
CALL  01 SPL  A (MSG  » 3 1 1 3 

readc i ♦*)  epsmin.rho.eta»kappa.oelta 

...DESIRE  TRACE  OPTION  f 

CALL  DTSPLAC9HTHACF  T/F,l,l3 
RE  AO  C i « 390  3 I TRACE 
FORMAT ( A 1 ) 

TTCITRACE  .NF,  1 HT  3 GO  TO  8? 

TRACE* , TRUE . 

GO  TO  83 

IF  C I TRACE  .NE,  IMF  3 GO  TO  81 
trace*. FALSE. 

,, .Input  functional  epsilon  ffps 
CALL  0 1 SPL  A ( 8HF  C T EpS*tl»l) 

PE  AD ( 1 * * 3 FEPS 

CALL  RAZOR (PH0*EP3MiN,PH0tFTA.KAPpA) 
wR I TE (6 . 39 1 ) (PHO(I) # 1*1 «K3 
format ( BE l 5 , 8 ) 

on  to  \ 

ENO 

FUNCTION  UCNAMF.PH03 

DIMENSION  PMO(l) 

common  /RZR/  K,0ELTAiALPHA.LlMlT»N«N?»N3«EPSL0N,THFTA(153f 

* 3(1*53*  FFPS  »TRAcF»LOW  (153  *HIGH(  IS  3 .RANGE  ( 1*5  3 
LOGICAL  TRACE.OOO.NqTIST 

REAL  LOW 

COMMON  /UFCT/  0(133, FPSno3,NS*NSl  ABS.1P.OOD.fi  O.FHI.FINCR, 

* OBI »SNl .SNINC1 »0P?fS^2*SNINC? 

comple*  rt 

Data  TWOPI/6.2831853071796/ 

...TRANSFFH  PERTURBED  VALUES  PpPS ( 1 3 1*1  TO  IP*  F RPM  PhO  VECTOR 
...TO  APPROPRIATE  EPS  VECTOR  ComPONFmTS  FOR  SUBROUTINE  RXTX 
DO  2?  T * 1 ♦ t P 
Tx*2*I 

FRS(IX)sPmO(I3 

T F ( ODD  .AND.  (I.EQ.tP})  GO  TO  ?? 

T x*2* ( I P*  I ) 

Iv*NS/2*l -I 
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f ps ( i x ) «»Hp( r v ) 

22  CONTINUE 

C transfer  pewtupbEO  values  t«IP*i  to  kb  fro*  PhO  vector 

C • • • • « T O D VECTOR  *I«1  To  NSLABS*\V  F-OR  SUBROUTINE  RXTX 
NNaNSu  ABS» 1 
00  25  T»i»NN 
T x * I ♦ I p 

23  DfI)*PHC(IX) 

NfR*InT(  (FhI«FlO)/fINCR*0.55  ♦ 1 

ITmP«inT(SN1/SNINCW0.5)  ♦ l 
NFCfcNT* (NFRtl ) /2 
OFPSET«10.0**L ,05*DB1  ) 

N0T1ST«. FALSE . 

00  26  I s l » ! ThP 
SKN*(l-n*3N!NCl 
00  27  J«1 *NFR 

TPltaTwOPI*(FLO*(J-l )*F1NCP) 

CALL  RXTx(NSLABS»SKNi»TPIL*CP8»0,BT) 

IJMNaC  ASS  C RT)-0PFS£T 
IFtNOTlST)  GO  TO  28 

notist».tpue, 

U^T  N*UMN 

26  IF  ( UMN  .LT.  IJMIN)  UmINbumn 

27  CONTINUE 

26  CONTINUE 

T Js-AmTNI (UMJN#0,0) 

- 

U* I Ns  1 0000  • 0 

0' - SET*10.Q**( , 05*0r2 ) 

^OTISTb.FaLSF, 

SkNbSn? 

100  V«UMlN 

00  2P  J*1 .NFP 

TpIl»TW0PI*(PL0t(J-1 ) I NCR ) 

CALL  RXTX(NSLABS»$KN.TPIL.EP3.0,PT) 

TFCJ  .FQ.NFCENT)  UMlN*CARS(PT) 

UMXaCABSfRTJ-OFFSET 
TF(NOTIST)  GO  Tp  30 
NOT  1 S T * , TRUE . 

Lima  x*uM* 

30  JFfU^X  , G T , UMAX)  UmAXkUMX 

29  CONTINUE 

Sk  NX  Sk  N*  SN I NC2 

IFIUMIN  ,LT,  V)  GO  TO  100 

T2*AMAX1 fUMAt,0,0) 

UbT 1 ♦ T? 

IF  (TRACE)  aRTTE(6*102)  NAME  ,IJ»0e  LT4  »FPSL0N«  (PhO  ( T j » T«1  *»0 
102  FORMA T Cl X.A6.UX,3G20.10/f 1 lx*5G20.10)) 

RETURN 

Eno 

subroutine  oispla (MSr,f Nwns. lines) 

OIMENSION  MSG ( 1 ) »MC g) 

OATA  NCRLE/  6001 OOOoOOOOOOOoOOOop  / 
IFL0(l*K,]lN)a3HIFTClw»I»K),AN0.SHlFT(MASK(K)tK) 
NxMAX0fMTN0(N*DS«6)  , 1 ) 

00  100  Isl.N 

ioo  Mfi)«MSG<n 

msuBNbM(N) 

00  102  lal.S 
J«12*(B»T ) 
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TFC  !FLD(J» 12.MSUBW) .NF.2R  ) GO  TO  jO« 

1 0^  CONTINUE 

10«  ! F f J ,LE.  36)  GO  TO  10? 

N*N*1 

*(N)«nCRLF 
GO  TO  106 

105  JbJ+12 

M(NJeuR(  ANO(MSUHN,MASK  (J)  1 ,SHIFT(njCPlF.-J)  ) 

106  NsNtl 
w f Kj ) so 

TF(lInFS  ,F.O,  2)  *BiTgf6»107) 

107  FORMAT ( 1 H ) 

WRITE ( 6 ) (M(T)«Ial.M) 

PpTuPNi 

END 

•6 
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PROGRAM  Tfc  ST?  ( INPUT,  TAPE.  1 ♦ OUTPUT  , T *Pf  fe»0UTPUT  * TAPE  1 0) 
C»»..,ThI6  PROGRAM  C^MPuTfS  THE  OIiARTfR  NAVE  LENGTH  DESIGN 

C TT  IS  USED  IN  CONJUNCTION  *iTh  R*S  ♦ RZ83 

COMMON  /R?R/  k .delta  . Ai  PHA  . I IVT  t ,NfN?»N3»EpSL0N,  theta  ns)  » 

• SU5)  .FEE'S.  TRACE  »lOI*f  15)  t HIGH (15)  .RANGE  MS) 
LOGICAL  TRACE. ODD 

Rf  -L  LO* 

COMMON  /DECT/  P ( 1 i)  .EPS  C 1 a)  ,N3iNSl.  AB9t  TPtOOO.El  O.FHI  .f  7NCR» 

* DB1 »Sn! ♦ 8 N I N C 1 .DR?. s^2, SN INC  2 
D!MfcNsiON  MSG(3) .MSG?(?) »hSg3(2)  . PhD (1 5 ) • 0 1 ElC T ( 1 0 ) 

DATA  mSG/1 Onf PSMIN.PHO  . 1 Oh.E T a «k APPA  « 7H.DEI  TA»  / 

DATA  wSG2/!0mFBEO  l 0 • h T • bH.INfR*/ 

input  ns.  the  number  oe  non. air  dielectrics 

.....NSLARS  IS  The  TOTAL  NUM8ER  OF  A)R  X NON-AIR  DIELECTRICS 
TP  IS  ThF  NUMBFR  op  DIPEERP.T  NON. AIR  DIELECTRICS 

is  the  dimension  oe  pho  vector  #totai  nr  oe  perturbations* 
call  DISPla(Rhn«  Si aRS* » 1 » 2 ) 

RE  AD ( 1 * * ) NS 
TPCNS  ,LE.  0)  STOP 
N9L ABS»2*NS 
NSLARi«NSlABS-1 
TP*(NStl)/2 

.... .DETERMINES  ip  ns  is  odd  or  even 

SETS  KsIP*NSi  A.01-2  IP  NS  IS  EVEN 

.....  * *eIP*NSLABl-l  IF  NS  IS  ODD 

.....thus  reducing  the  number  op  perturbation  components 
.....on  The  Pwp  VECTOR  As  a result  of  the  MIOOLE  OIElFCTRIC 

VARYING  AT  l.AMBOA/(«*SQRT(fP$IlON)) 

nooa, TRUE, 

K*IP*NSlABi •) 

IF (mOO(NS.2)  ,EU.  1)  GO  TO  a? 

COD* * r A u SE • 

K B X • ! 


C TNPU*  DIELCT(I)  1=1  tq  IP  t JNTTTaLI7E  EPS  VECTOR 

ENCODE  (20.28T.msG3)  IP 
26R  FORMAT!*  EPS(1  TO  *,T2«*)«*) 


CALL  DISPLA (MSG3.2. 1 ) 

RE  AQ ( 1 . * ) (OTELCT ( I ) . 1*1 . IP) 
DO  7P  1*1 , IP 
I x*2*I 

EPS( T X • I )*1 ,C 
F PS ( I X ) *0] ElC T ( 1 ) 


IE (ODD  .AND,  (I.EU.TP))  GO  TO  7* 

T x*?* ( TP* j ) 

F»SU*»1)*1  .0 

I y *nS/2«  t -T 

EPS( I x)*OIF  LCT ( I V ) 

T«j  CO*  T TMIfc 

C Input  0(i)  i*i  t<  nslari 

ENCODE (?r ,290. M5G3)  NSLARl 
2R0  EoRma 7 t » D ( 1 TO  * » T 2 . * ) ** ) 

CALL  DTSPlA (MSG3.2. 1 ) 

READ ( 1 .*)  (0(1). I»1 ,NSl  APj ) 

C INPUT  CONSTRAINTS  l OH. HIGH. RANGE  ON  ENTIRE  PHfj  vFCT()P 

C INPUT  LOh(T)  l*l  TO  k 

encode  (2o  ,?<n .msgs)  k 
29l  FORMAT!*  lO*(t  TO  *,I2.*)r*) 

CALL  OlSPl A (MSQ3.2.  i ) 

RE  AD ( 1 * * ) (LOrt(I) tT*l  .K) 

C Input  HlGH(I)  T*1  TO  k 
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ENCODE(20.292,mSG5)  K 
292  F0«M4T(*HIGH( 1 TO  F.T2.*)**) 

CALL  Dl3PLA(MSG3,2i i ) 

RFAOfi,*)  (HIGH ( I J *i«i  ,K) 

DO  8A  I ■ 1 * K 

8«  RANGE ( I) ■100,0* (HIGHf I)*LOWf I) ) 

c INITIALIZE  PHO  VECTOR  to  INPUT  DIELECTRICS  X THICKNESSES 

DO  79  1*1 . TP 

79  PmO(I)*DIELCT(!) 

• NN1«N3LAB1/2 

NN2*NN1 42 
NNJsNNt 

! J*  1 

DO  80  1*1 ,NSLAfll 

IP!J»IP*IJ 

IF ( ODD ) GO  TO  801 

IFd.EQ.NNl  ,OR.  I,rO.*N2)  GO  TO  *0 
PHO(IPTJ)«D(n 
T J»I Jtl 
GO  TO  80 

801  I F ( I ,EQ.  NN3)  GO  To  80 

PHO(rpijj*Dcn 
T J*  I J ♦ 1 

80  CONTINUE 

C Input  frfqloh.freqhi .freoincR 

CALL  DTSPLA (M3G2f2. 1 ) 

R£ADd»*)  FLOtFHl »FlNCB 

c Input  obi 

CALL  DISPLA(«HDbt*M»n 
PfADfl**)  OBI 

C. ....INPUT  SIN(THETAI)  X SIN(ThFTA1 JlNcRfcMENT 
CALL  OISPLA(10HSTHltlNCR*»1  .1) 

RE  AD  (1  * * I SNl.SMNCl 

C Input  OH?..,,,, 

CALL  DTSPLA  (<»mDH2«*  1 • 1) 

RFADd«*)  D«2 

c... ..Input  SIn{Thfta2)  X SIN(ThFtA2)!NcREmENT 
CALL  DISPLA(10HSTH2XINCR*»! ,1) 

Rf  AD  f l * * 1 SN2.SNINC2 

c Input  alpha  x max  lihit  of  u function  evaluations 

Call  DTSPLA(9HALFA.LwT*.l .1 ) 

BE  AD  ( 1 t * ) ALPHA. L I M I T 

c. ....  INPUT  F PSM I N.RwO.ETa. kappa. DELTA 
call  dtspla (msg, 3, n 

Pf  ADf  1.*)  EP3MlN»RHf).ETA.KApPA.PFl  TA 
C..,., DESIRE  TRACE  OPTION  f 

81  CALL  OTSPL A (OhTRACF  T/F.I.tl 
RE  AD ( 1 . 390  ) T TRACE 

390  FORMATfAl) 

IF(ITRACF  ,NF.  1 H T ) GO  TO  ftp 
TRACE*. TRUE, 

GO  TO  83 

82  IFdTRACF  .NE.  1 HF ) GU  T.O  *1 
TRACE*. FALSE, 

C.....TNPUT  FUNCTIONAL  EPSILON  FFPS 
es  CALL  DISPLA(8HFCT  EPS«.1.11 

RE  AD  C 1 . * ) FfcPS 

CALL  RAZOR fPHO.EPSMlN.PHO.FTA .KAPpA  ) 

WRITE  ( 10. 391)  ( RHO  d)*T*l.K) 

391  FQRha  T ( 5F 1 5 , 8 ) 
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GO  TO  1 
E NO 

JUNCTION  U(NAME.PHO) 

DIMENSION  PHO( 1 ) 

COMMON  /RZ»/  KtDELTAtAl  PHAtLTMlT,N*N2iN5.EP8LON,THETA(15)  » 

* S(15)»FEPStTRAcE*L0W(15)»HIGH(l5)*RANGP(l5) 
LOGICAL  TPACEtOOOtNoTlST 

PfAL  LOw 

COMMON  /UFCT/  D(1 J5 ,EPS( 1 a) ,NS»MSl ABSiIPf ODO.FLOtFHl »FINCR. 

* OBI »SNl »8NINC1 »DB?,SN2tSNlNC2 
COMPLE*  BT 

DATA  TWOPI/&.2831B55071 79b/ 

C TRANSFER  PERTURBED  VALUES  »EPS(n  I«1  TO  IP*  FROM  PhO  VECTOR 

C,,,,.To  APPROPRIATE  EPS  vFCTOP  COMPONENTS  FOR  SUBROUTINE  RXTX 
00  22  IB1 *IP 
TX«2»I 

EPS ( I X ) *PhO ( I ) 

IF  ( ODD  .AND,  (I.EQ.iP)I  GO  TO  22 

ix«2*cip*i) 

IVaNS/2+I.I 
EPS ( I X 5 sPHO ( I Y 5 


22  CONTINUE 

C,,., .TRANSFER  PERTURBED  VALUES  *0(15  IaIP*l  TO  <•  FROM  PhO  VECTOR 
C*..,.TO  D VECTOR  #1*1  TO  N8LABS-1*  FOR  SUBROUTINE  RXTX 


PwniPBPHO(iPl 
NNaNSl ABS»1 
nnI «nn/2 

NN?"NN 1 
N K 3 a N N J ♦ I 
T J«i 

00  2 3 I ■ 1 » N N 
IJTP«IJ*TP 
TFfQOoI  GO  TO  20 

IFC.E0.NN1  .OR.  I.E0.NN2)  GO  To  \<t 
Of! )BPMOf I JIP ) 

T J*3 J*! 

GO  TO  ?3 

IV  0(1 )«0.25/SQRT(PHOTp! 

GO  TO  23 

20  I F ( I , F 0 , NN3)  GO  T0  IP 
Of  T )apHQ( I JIP) 

IJ«UM 

23  CONTINUE 

nfr* I NT  ( (FHI-FLO)/FTNCR+O.S)  * 1 
ITHR«INT(SN1/SNINCU0,5)  ♦ 1 
NFCENTa(NFBt1 3/2 
nFFSET*10,C**(.OS*nRi  ) 

NOTlSTa. FALSE. 

00  2ft  I * 1 » I T MR 
S*  N*  ( 1-1  3*8MNC  1 
00  27  Ja 1 .NFR 

TPIL«TWOPI*CFLO+ ( J-i )>r I NCR) 

CALL  RXTx(N3LABS.SKN.TPIL.EPSfD,RT) 
UmNsCaBSCRTJ.OFFSET 
TF(NOTIST)  GO  TO  2F 
NOTlSTa. TRUE, 

L)  m I \ a i j M n 

26  IF(uMn  ,L  T,  IIMIN)  nvINauMN 

27  CONTINUE 

26  CONTINUE 
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!JmTN*1  0000 ,0 
OFF3ET»10.0**(.0S*Db2) 

NOTlST*. FALSE. 

9KNaSN2 
100  V«IJMIN 

no  29  Jai.NFR 

TPIt8TwOPI*(FLn+f J-i )*FJNCR) 

• CALL  RXTx(N3l ARStSKN.TPIt.EpS.O.RT) 

T F ( J .FGl.NFCENT)  UM I N«C ABS ( RT ) 

LJMX«CAR8(RT)-0FFSET 
TF (NOT  1 ST ) CO  T(j  30 
NOT 1 ST» , TRUE . 

UMAX»UMX 

30  I F ( X ,GT.  UMAX)  UMAXBUMX 

29  CONTINUE 

SkNsSkn*9NTNC2 

IFCOMIN  ,tT.  V)  GO  TO  100 
T2*AmaxI (UPAXtO.O) 

U»T  t ♦ T ? 

TFtTRACE'  *R I TF  ( t>  1 1 02 ) N AME  ,1),  DFL  T A #£P8lON,  (PHO  ( I ) . I « 1 . K) 
102  FORMAT (lX.A6«aX,3G20.10/( 11 X»5G20. 10)) 
return 

F NO 

subroutine  displa (MsG*Nkvos.LiNFs) 

0 T Mfc  N9 1 ON  MSG(1)»M(«) 

Data  NC«IF/  OOOIOOOOOOOOOOOOOOOOB  / 
tFlD(ItK»IW)«SNlFT(Iw,I*K).AND,SHlFT(MA8*<(K)«K) 
N«MAX0(MTN0(N^DSi6) ,1 ) 
rO  100  I ■ 1 • N 
100  -(I)aMSG(I) 
msubn»m(^) 

DO  102  Ial.S 
Js 1 2* ( 5*  T ) 

I F ( TFLO(J.  12,M8UtjNi  .NF.2R  ) GO  TO  10<J 
102  C DN 7 * Nllfc 

1 0 u IFCJ  ,lE.  3b)  GO  TO  105 

N«N  ♦ 1 

m ( N ) snC^LF 
GO  TO  1 ON 
105  J*J*12 

Mf*.;*(jRUNC(vSwbN»w6  3K(J))»$MlFT(NCRlF**J)) 

10b  N ■ N ♦ 1 

M(N)»0 

I F ( l I NE  S ,FQ,  2 j *RlTE(b.l07) 

107  F “flv  A T ( 1 « ) 

*R I T F ( b ) (M(l)*I»l*NJ 

rfTupn 

END 

• • 
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n n n r»  noon 


PROGRam  Tf  5T?(  INPUT  .TAPE  1 ♦ OUTPUT  . T APE6»0UTPUT , TaPE.  1 0 ) 
.....This  prog  am  IS  ESSENTIALLY  TEST1  EXCEPT  that  IT 

.....incorporates  a ta&le  search  for  physically  realizable 

.....OIElECTRIC  MAYfRlAUS.  IT  IS  USED  IN  CONJUNCTION  HlTH 

r x s ♦ RzsafrtHicn  nots  a table  search  in  subr.  fxplR) 

COMMON  /RZR/  K.OElTa.AL PHA, limit, N.N?»N3,EP8L0N, ThETA(15) . 

• Sf 15) » F EPS. TRACE. LOW ( 15) ,HICH( 1 5 ) .RANGE ( 15) 
LOGICAL  trace, ODD 

real  lOw 

Common  /UE'T/  D(13),EPS(lR),N3.NSt ABS»IP.ODO,FLO,EHI.EINCR» 

• OBI  *SNl  ,8MNC1  ,DB?.SN2.SNINC2 
DIMENSION  mSR(3) .M8G2C2) ,MSG3C2) . PHO ( 1 5 ) , 0 1 ElCT ( 1 0 ) 

DATA  MSG/lOrfPSMjN.RHO  , IOh.ETa, kappa  , 7H , DEL T Ac  / 

Data  m$G2/ i Ohfreq  lo.ht  , INcr*/ 

input  ns,  The  NOMafp  of  NON. AIR  DiELFCTRICS 

NSLAPS  IS  The  TOTAL  NUMBER  OF  ATR  X NON. AIR  DIELECTRICS 

IP  IS  THE  NUMBER  CF  OIFFFRENT  NON. AIR  DIELECTRICS 

. . . . . K IS  ThE  DIMENSION  qE  PhO  VfCTOR  #TOTAL  NR  OF  PERTURBATIONS* 
1 CALL  DISPlA(RHNP  SLABS*, 1,?5 

READCl,*)  NS 
TF(NS  , LF » 0)  STOP 
NSLABS*2»NS 
NSL  ABl«NSLABS-l 
jp*(Ns+n/2 
K«IP*NSLABi 

c Input  DIElCT  f I ) I«!  TO  ip  X initialize  EPS  VECTOR 

FnC0DEC20.2BO.msG3)  TP 

289  FORMATf*  EPS ( 1 TO  *,T2.*)b*) 

Cii.L  DI3PLA  (MSG 3. 2,  { ) 

PEADC 1 ,*)  (DIFLC  T ( I ) » I«i , IP) 

000*. TRUE . 

IF  fMOO(NS,2)  ,ECi.  0)  OODB.FftLSE. 

DP  78  Isi  , TP 
T x«2*  I 

FPS{IX-1  )el  ,r»2 
FPS(IX)=OIFLCT(I  ) 

IFfODD  ,ANO,  (I.MJ.fP))  GO  TO  78 
T x«2*  HP*  I > 

FPS( J*-i ) s i .02 

I YB\S/2» 1 -I 

FPS( IxJsDIFlC TCI  Y) 

78  CONTINUE 

C INPUT  D(T1  I«1  Tf)  NS!  API 

FNCCntf20,?RO,MSG3)  NSLAPj 

290  FORMATf*  DM  TO  *,T2.*)«*} 

CALL  OISPlA  CMSG3-2. 1 ) 

READfl,*)  fD(I).I*l  ,NSL  ABM 

C INPUT  CONSTRAINTS  1 r*1, MIG*', RANGE  nN  ENTIRE  PHO  vFCTOR 

C.. ...INPUT  LOh(I)  lei  T0  * 
encode (2o.?Ri ,m«g3)  * 

291  rnevATf*  L('»li  T 0 *,T2.*)b») 

CALL  DTSPla (“SG3.2.  \ ) 

Rf  AD(  1 * * ) f LOwC  I ) ♦ Tsl  .tO 

C TNPoT  Hlr.nm  1*1  To  K 

FNCuOE (20.?9?.msG3)  * 

292  FORMAT  (*hIGh(1  TO  *,T2, »•)«*) 

CALL  DTSPl A (MSG3.2.  | ) 

RF ADf i .* ) (hIGh( I ) , 1*1  .k  ) 

DO  BA  T*1 

8<*  range  (T)bi oo,o*(Hic.Hf  n-Lo**fi) ) 


C - i (J 


C INITIALIZE  PHO  VFCTnfi  TO  INpUT  OlFuECl»ICS  X THICKNESSES 

PO  79  1st  , TP 
79  PhD ( I ) »D TfcLC  T ( I ) 

00  80  Ib1,N3lAB1 
TX«IPtI 

sc  PHO(ix)ap(n 

C.,,. .INPUT  FPFOLOW.PPEQMI ,f RE0INCP 
CALL  D T SPL  4 ( WSG2  • 2 « 1 ) 

»EAl>d»*>  PLD.FHI ,FINCP 

C,,,,,lMpuT  081 

CALL  OISPL A (UHPH1«.  l .1) 

PEA0(1«*I  081 

C INPUT  SINCTHFTAO  X SINCThETAniNcPEMENT 

CALL  OISPLAClOMSTHlxTNf«»,l #l ) 

SN 1 » SN  T NC l 

C., ,,.TNpuT  Oft? ... 

CALL  OISPL A (aH082»* 1 » n 
PEADd»*>  OB? 

C INPUT  9IN(ThFTA?I  X SIN (ThFTA?> increment 

call  DlSPLAUOHSTH?x!NCR».t  ,1) 

RFA0(1.*T  SNP.SNINC? 

C input  alpha  X MAX  L I M I T OF  u FUNCTION  EVALUATIONS 

CALL  DTSPLA(RHALFA,LMT«.l ,1 } 

BEADCl.*!  AlPHA.LIMjT 
C. .. ..INPUT  FPS^IN.RHO.ET a. KAPPA. DELTA 
CALL  OISPL A (MSG. 3. I ) 

RFADCl  . *)  EPSMlN.RHo.ETA.KApPA.DFl TA 
C.,.,,PESIRfc  TRACE  OPTION  f 

81  CALL  DI SPl  A (RhTRACF  T/F.1,1) 

PEAOfl.390)  ITRACE 

390  FOPMAT(Al) 

TFflTRACF  , N t , 1 H T ) GO  TO  8? 

tracf«.truf, 

GO  TO  81 

82  TFflTRACF  , Nr  , JHF)  CiO  TO  8\ 

TRACE e. FALSE. 

C ..... T UpuT  FUMCTIONAI  E P S 1 1.  P N FFPS 

83  CALL  DTSPLAfSHFCT  EpSb.1.1) 

R F A 0 ( 1 . * ) FEpS 

CALL  RA?PRfP80»EPSMjN»RH0»FTA.KAPpA) 

*«T  TF  £10. 391  • f PhQ  ' i ) , l«i  ,kj 

391  FORMA  T lS.rl) 


r.n  TO  l 
F NO 

function  i j ( n a m e » p h n i 
dimension  PmO(1  ) 

common  /RZR/  K.OELTa.AI PHA.I  TMIT,M.NP#N3«EPSL0N,THFTA(151 . 

• Sf  IS)  .FFPS.tr  ace  .LPW(1S).HIGH(  IS)  .RANGE  (1*5) 
LOfWCAl  TRACE  .000.  NoTlST 

REAL  LPft 

roMMON  /OF  C T / 0(1  i).FPS(U)  .NS.nSl  AHS*TH.ODD.FU),Fhi  .FTNfR. 

* 0H1 «SNl .SNlNCl .DR?.SN2,SNINC? 

COMPLEX  R T 

0 a T A T*OPI  /fc.?8518\>3071  79b/ 

C.....TRANSFFW  PERTURBED  vH"ES  PFPSd>  Tel  TO  IP*  FROM  PhD  VETTOR 
0 .....  T 0 APPROPRIATE  f PS  VECTOR  ComPOKiFnIS  for  SUBROUTINE  rxtx 

on  d?  1*1  .ip 

Tx=<?* ; 

fPSC  Tx  HPHH(  T ) 

IF f ODD  , ft NO,  Cl. EG. IP))  GO  TO  ?? 


C-ll 


TX«2*( IP*I) 

I y«N8/2* 1 • I 
E'°s  ( I x ) »pho  ( t y ) 

22  CONTINUE 

c . . . . .TRANSFER  PERTuRRED  VALUES  *0(1)  T k T P4- 1 TO  KF  F«0M  PhO  VECTOR 
C.....TO  D VECTOR  #I«1  TO  NSLA8S*i*  F(iB  SUBROUTINE  RXTX 
NN«NSlA0S»1 
00  2 3 T»1 * NN 
T X*  I ♦ I P 

23  0 { I ) *PhO ( I X ) 

NF  R* I N T ( (FmI«FlO)/fTNCR*0.S)  ♦ 1 
!THRaINT(SM/SMNCUO.S)  ♦ \ 

NFC£NT*(NFR+1 )/2 
OEESET«10.0**(,05*ORt ) 

NflT  1 9T» . F AlSF  . 

00  26  1*1 » T TmR 
SkN«CI«1  )*SMNC1 
00  2T  Jm\ ,NFR 

TPIl«T^OPI*(Fl0«(J-1 )*FINCR) 

Call  »XTX(NSLABS*SXN»TPIL tFPS*D,RT) 

UmN*CaR3(BT)*0FFSET 
IF(NOTIST)  GO  TO  26 
NpT 1 ST« , TRUE . 

UM  I 

28  IF(UmN  .Lt.  UMINJ  U^INbUMN 

27  CONTINUE 

26  CONTINUE 

?1 a-AMlNt (UWIN»0,0) 

IJMI  N*  1 0000 . 0 
r)FPSFTal0.Py*(,05*r>R2) 

NOTlST*. FALSE. 

3*  NsSN2 
100  Vs'jMlN 

DO  2R  J s 1 * NF R 

TPlL«TwOPI*(FLO+(J*l)*FINCR) 

Call  rxtx(nslabs*skniTpil*f pSiO.pt) 

IF(J  .EO.NFCFNT)  U^jN.CAPSfRT) 

UMX»CA«SfPT)»OFFStT 
I F f NOT  1 ST ) GO  TO  3n 
NOTlSTc.TRUf . 

UM A X *UMX 

30  IP(UMX  ,GT.  umAX)  UmAxbuMx 

29  CONTINUE 
Sxns8kn»SnInC2 

If f um i n ,lT.  V)  GO  TO  100 
T2*AMA  X 1 (UMAX«0,0) 

U«T 1 ♦!? 

iFfTRACE)  ^RITE  (f>.  1 02)  N A*<F  , |j  , OF  L T A i EPSLON  * C PHO  ( I ) « I ■ 1 . K ) 

102  FORMAT ( l X . Ab.PX ,JG2o, 1 0/(11 x«S&20. 10)3 
B J T uR \ 

FNO 

subroutine  dtspla (msg*n*»os«l  tnesi 
OI-ENSION  M8r,(i)*M(EO 

oata  nc^le/  ooo i oooooooooooooooop  / 

TFLOU  *«  * IX0*3MIET  f iw,  T*K)  , aNO.SMIET  fMASK  (K) 
N*MAXO(MINO(klx(OS«6)  . 1 ) 

00  100  I a J » N 
100  M ( T ) *mSG ( I ) 

M3UoN»M(N) 
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np  i <v  i«i »s 

jai?*cs-n 

I F ( IF  LQ ( J * , NF  « ?W  ) GO  TO  1 OU 

10?  continue 

100  I F t J .It,  561  GO  TO  10*5 
N»Nf  1 

M(N)sNCRLf 
GO  TO  106 

105  J«  J ♦ 1 ? 

( 4NO(M3URN,MASK ( J) } .SHIFT  ImCRLF * • J ) ) 

106  N»Ntl 

M(N)S0 

IF ( L I NP S , F Q , ?)  aPiTE(6.107T 

107  FCRMATdH  ) 

WP  T TF ( 6 ) (*(Y)«I»1«N) 

return 

ENO 

M 
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C.,..,ThIS  SUB-PROGRaM  IS  CALLED  R*S 

C.....TT  IS  USc 0 IN  CONJUNCTION  *iTH  FlTHER  TE 3 T 1 / T E ST2 / TE S T3/ 
SUBROUTINE  RXTX(NSLaH*9KN» TpIL»EPS*0»RT) 

COMPLEX  RT*Ut  J*GL.RT,TY,8(2*2)  *A(2,2)  ,AD(1«*2#2) 

DIMENSION  FP$(l)*D(l),Y0(14).AJ(U.2,2) 

RFAL  kAPPA(IO) 

FUUIVALENCE  (SlliS(iil))  , (321,3(2*1))  • (RLTY,?Y) 

DATA  J / ( 0 , « 1 . ) / 

c 

COMPUTt  kaPPavS  and  YC VS  #PFAL* 

c 

3kN2»SKN*5XN 
DO  100  NOsi.NSLAS 

« APPA  ( VO)  sSQR’’ (EPS  (KjD) -SKN2  5 
100  YO(ND)*EPS(ND)/KAPPa(ND) 

c 

C CONSTRUCT  JUNCTION  wave  TRANSMISSION  MATRICES,  AJ,  *REAL* 

C 

NSL  AB  l aNflL  Ab*»l 

DO  110  Np  b 1 , NSL  A B 1 

BLTYsI  ,0/(Y0(ND)*’Y0(‘ND*l  ) ) 

311*(Y0(ND)»V0(ND*1))*RlTy 
321«2,C*Y0(ND)*RLTy 
RlTybyO (ND  + 1 ) / YO ( NO ) 

A J ( ND  , 1 *1 ;al  ,0/321 
AJ(ND,  1 , 20*31 1/S21 
AJ(ND,2, 1 )8AJ(ND* I *2) 
aw’(N0,2.?)*(BlTY»S2i*S21*S1  1*31  1 )/82l 
C 

C CONS'WUCT  LINE  LENGTH  WAVE  TRANSMISSION  MATRICES*  AO,  BCOmPlEX* 

C 

AOCND,  1 , t )-*CExP(  JAKAPPA  (NDm  )*TPtL  *0(ND) ) 
ad  (N'1, 2<  ?)3  i ,/AD(ND,  1 , 1 ) 

AD(N0,1,2)s{0?*0.) 

110  AD(ND,2,1)*(0,,0.) 

C 

C FORM  TOTAL  *■  A VE  TRANSMISSION  MATPlX,  A,  BCOMPlFXB 

c 

DO  1 30  noc i , 2 
DO  i 3 0 mob',,? 

A ( NO  * MO ) s ( 0 , ,0.) 

IF  (ND.fcG.MD)  A{NJ,mO>b(1.,0.) 

130  CONTINUE 
C 

r o 1*0  ND* 1 * NSL  aB l 
ND1 *Nsl  AR-ND 
DO  ISO  I Os l ,2 
DO  IbO  K0sl,2 
T v * ( 0 , ,0.  ) 

DO  1 5 C L D* 1 ,2 

U* (0, , 0 , ) 

00  1«0  m p a 1 ,2 

140  IJbU*AD(ND1  ,lO*mD)*A(MO,kD) 

ISO  TV*TY*U*Aj(NDl,in,LD) 

100  S(IO*»<n)sTY 

DO  170  I D 3 1 ,2 

no  no  KO«l  ,2 
170  AtID**tn)B8(TD,KD) 

ISO  CONTINUE 

C,tHTB(YO(NSLAH)»YO(j))/(YO(NSLAR)»YO(l)) 
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T Y « A ( 1 « n ♦f.L  B T * - 1 1 

HFTuPNvS  TRANSMISSION, 

RT»C1  .♦&C«’r  )/TV 

retoRn 

emd 


• F tc  lf*t 


#CDMPLt*» 


C-15 


C.....THIS  SUBPROGRAM  IS  CALLFD  RZS3 

C.....IT  IS  U8ED  IN  CONJUNCTION  with  EITHER  TFST1/TEST2/ 
SuRPOuTINE  RAZ0»(PM0tEPSMINtRH0,ET4*KAPPA) 

DIMENSION  PMI(15)»PH0(1S) 

COMMON  /BZ»/  K.DELTa*  ALPHA,  l IMTT  » N » N? . N3 » EP8L0N , ThE'i  A (1  5 > . 
* S(15) *FEPSt TRACE* LOW (151 .HIGH  05) .RANGE (15) 

LOGICAL  FINISH. TRACE 
REAL  L'Ow 

call  sECOnD(TIme) 

Ns  i 
N 1 * l 
N2S0 

njbO 
NusO 
I C * o 


C Evaluate  function  at  starting  poTnT 

UPM0«U(6HRAZ0R1 |PHO) 

IP(TRACE)  wRlTEffe.l)  U P H 0 • ( P H 0 ( T ) • I ■ 1 *K) 

1 FORMAT(lHO.ay.6HUPMo  *.G20.10«15*,5HPHO  ■ ♦ 3G20 . 1 0/ ( 1 1 X ♦ 5G20  , i 0 ) I 

C FINISH  CRITERIA  SATISFIED  f 

IF  (FINISH(UPHO)  j GO  TO  1 <1 
DO  2 I * 1 |K 

2 9 ( T ) * 1 , 0 
EPSLONsEPSMIN*  CETA**KAPPA) 

C,.,. .CONDUCT  PATTFRN  SEArCH 
IFCTRaCE)  WRITE(ft.loO) 

100  F0»MAT(7H  PAT3R1) 

CALL  PATSER(PHOtUPHO) 

C FINISH  CRITERIA  SATISFIED  f 


IF (FINISH ( UPHO 1 ) GO  TO  l <i 

00  10  J11, KAPPA 

NaaN»N41 

UPH0«U(6hRAZ0R2«Ph0) 

IF ( TRACE  I wR I Te (6.31  N4»UPH0» (PHp(I) » I ■ l * K ) 

J fhRmA T ( Uh  n s. 15. 16X.6HUPH0  ■ . G20 . 1 0 /6X . 5HPH0  S.3G20.10/ 

* (11X.5G20.10)) 

NasN 

!PfN.cT,LIMIT)  GO  To  1? 

c.,..,ortain  new  starting  point  and  exploratory  parameters 

OElTasO  ,0 
DO  u I*  l ,K 

RanDOms-i  ,0f2,0*RAMF(0,0) 

Pm  I ( I )«B0UND( I *PHO( I )#RH0*RANOOM*FPSL0N/100t0) 

TstPHj  m-PHo(l) ) /range (!)*i oo.o 
a * oeltasdeltaatat 

del  TA.SGRT (DFLTA/k ) 

c . . . , ,f valuate  function  ano  conduct  pattern  search 

UPHl»u(6H«AZ0R3.PHl ) 

NsN  ♦ 1 
N 1 «N 1 ♦ 1 

IF ( TRACE  ) ARlTE(b.5)  uPhI  , (PhT ( I)  ♦ Isl  .K  ) 

FOPMATf IHQ.ax.bHUPMi  *»G?0. 1 0. 15X.5HPHI  S.3G20.10/ 

• (11X.5G20.10)) 

FPSlOf  sEPSLON/ETA 

TF(TRaCE)  wRITE(6.101) 

. ; Fr«jMAT(7H  P A T S R 2 ) 

'»  . PA TSE R (Ph! .UPhi ) 

....'■  AN  IMPROVEMENT  i 

(u»Nl .LT.UPhO)  GO  TO  7 

* .uT.lI“It)  go  to  i? 


016 


C.  ...  .PROJECT  ME*  POINT 
OELTAsO.0 
DD  6 1 * 1 * K 


T«PHO<I)-PmKX) 

TIbT/HANCF (15*100.0 
ThETA(I)sT 

PhI ( T )«BOUND( I .PHO(l )♦!) 
b DELTA* OElTA.T1*T1 

UPH I«uPhO 

C .... .CONDUCT  PATTERN  WOVE 

IFfTRACE)  WRITE (fe* 1 0?) 

102  FORMA  T ( 7m  PATmVI) 

CALL  PATWV (LPhI *PHI .PhD) 

C.,,, .OUTCOME  an  IMPROVEMENT  ( 
IECuPmI.GE.uPhC)  GO  TO  9 

C RETAIN  best  POINT  AMD  FUNCTION  VA|_UE 

7 UPH0«UPHI 

C . . . , .PROJFCT  NE*  POINT  o 


DELTAsO, 

DO  8 I ■ 1 |K 
TsPhI ( I ) -PMO  ' T ) 

T l sT /RANGE ( I) *100.0 
THE  T A ( I ) • T 
PhO(I ) sPh I (I  ) 

PhI ( I ) sBOUNDf I .PHO( I ) ♦ ! ) 

8 DEL TAsDELTA*T1#T1 

IF (N.GT, LIMIT)  GO  To  1? 

Cm,  .CONOUCT  PATTERN  move 

IFfTRACE)  wRlTECfetlQl) 

103  FCBMAT(7H  PATMV2) 

C-LL  PATmv (UPhI *PHT ,PhO) 

C .... .OUTCOme  AN  IMPROVEMENT  ( 

IF (UPhI ,L T .UPHO)  GD  TO  7 

c,.,, .finish  criteria  satisfied  f 

9 IF CF InISh(UPhO) ) GO  TO  1 R 

10  CONTINUE 

CALL  SECONDfTlMEi  ) 

T I ME  b T I Mp  i . T i Mf 

WRITE  (Fill  1)  Nf  TIMt 

11  FORMATf*ONO  CONCLUSION  BY  RaZO»  AfTER*. IU»*  FUNCTION  EVALUATIONS** 

• * A \0*  » F 7 ,2»*SECONDg*) 

GO  TO  16 

12  CaLL  SE  C OND ( T I me  1 ) 

T IME*T IMF l «T IMf 

wR I TE ( 6 * 1 3 ) TIME 

13  FORMAT  f5<ihOLlvlT  ON  FUNCTION  EVALUATIONS  EXCEEDED  B Y RAZOR  AFTER. 

• F7.2.8M  SECONDS) 

GD  TO  16 

U CALL  SECONDfTlMEi ) 

TIME*TTMF1«TIME 
wr I Tf ( 6 * 1 5 ) T l me  * N 

IB  FORMAT  ( 33HOCONVERGFMCF  OPT  A I N F.  D Ry  RAZOR  IN.F7.2* 

• l<iH  SECONDS  aFTER.IS.21  h FUNCTION  EVALUATIONS) 

16  NuaN»N« 

WRITE  (6.17)  N1  .N2.N3»N«»tJPHO.  (PhO(  I ) . I«l  .K) 

17  FORMATf37H  The  FUNCTION  EVALUATION  BREAKDOWN  1 S ✓ 1 1 0 . 

• 9h  BY  RAZOR/i  l 0 . 1 <IH  pY  PaTmV.  ANO/IIO.Rh  fly  E * PL  R / 

• 2UHCFINAL  PARAMETER  VALUES.*^  N a,I5»3x» 

• 6HjPhO  s.GlS.7,SK.SHpH0  «.g15.7/( 1 3X.6G15.7)  ) 

R E T ij  R N 
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END 

LOGICAL  FUNCTION  FImISw(UPHO) 

COwmOn  /R'B/  K, DElT4, ALPHA. lI^TT. n,N?*N3. £PSLON, ThETA(IS) * 

* S(  15)  , PEPS »r»ACE, LOW ( 15)  »HIGHU5),BAMGF(15) 

LOGICAL  TRACE 

<?£  A L LO« 

pr  1 sh*.  false, 

TE(UPhO.LT.FFPS)  FInIShx.TRijE  . 

IF(TBaCE)  **ITE(6,100) 

100  FOPvAT  f 7h  FIMSh) 

P F T U R N 

£ Kip 

Subroutine  patsebcpho.uphoi 

OIMfcNSlOM  PHI (15) tpHOf 1 ' 

COMMON  /»ZR/  K,OELTa, ALPHA, L TWIT, m,N?,N3,EP8LOM, theta ( 151 , 

* sn5)*FEPS,TRAcE*LPW  (151**1  GH  (15). RANGE  (15) 
LOGICAL  PI^ISh,  TRACE 

REAL  lOa 

1 T P C N , G T , LIMIT)  RFTURN 

c Exploratory  increment  too  small  r 

TP (DELTA, lT.EPSlON)  RETURN 
IE  (E  IMSh(UPhO)  ) RETURN 
UPWI *uPHO 

on  <t  i ■ 1 , k 

2 PhI ( I )*PHpri  ) 

C CONDUCT  EXPLORATION  around  pASE  PpINT 

CALL  £ X Pi.  R ( UPh  I , PH  I ) 

IP(T»ACE)  WRITECb.lCO) 

100  P0eMAT(7H  PATSER) 


C oltcomf  an  improvement  l 

rpfuPHl.LT.UPHO)  GO  to  3 

C REOuft  EXPLORATORY  jNCREMFNT 

7 DEL T AsALPhAADELTA 

GO  TO  1 


C...,.P£TAlN  BEST  POINT  AnO  EjNCTjON  V A | Ufc 

C project  nea  point 

3 DFLTAsO.O 

UPH(jaupH  T 

On  r I b 1 . K 
T 1 =PH  I f I ) -PhO  ( I ) 

TiT^/baNGF n 1*100.0 

ThETA(I)bT1 

PhO ( T ) BPHI  f I 1 

P-I1I) snOuND  f I , PHO  f T 1 ♦ T 1 1 


u DFLTiaPEl TA4T*T 

I F f F I N I Sh ( UPhO) ) RFTURN 
TF f N.GT.L IMIT)  RETURN 

C. ... .CONDUCT  PATTERN  move 

Call  patmv (uPhI ,rht .phpi 
TFfT»ACE)  w»TTE(6,100) 
c .... .outcome  an  improvement  t 

TE  f UpH  I »l  /PHO  1 3*7,7 
ENO 

function  bound u *e*pR) 

common  /PZR/  k,  delta,  alpha,  I TMTT,N,N?,N3,EPSLON,THFTAf 151, 


* Sf 15), EFPS, TRACE, LOW (151 , H I fin C 1 5 ) , R AnG£ (15) 

REAL  LOa 

50UNnsAMAxj(AMTM(Hir,H(I),ExP«),l0w(T)) 

RETljRn 

EnP 
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SuBHOuTlNfe  £*PLfi CuPhO.PnI ) 

DI^fcNStOM  PHI(i) 

COmmOm  /»ZR/  K. DELTa. ALPHA. lTMIT.N»N2. NS. EP8L0N,THETAn5). 

• S( 15) .FEPS»TRAcE.L0*(15) »HIGH( 15) .RANGE(lS) 
LOGICAL  trace 

Bf Al  Lnw 

no  « i*i .* 

SOEL*OELTA/iOO.O*S(l )*RANGE(I) 

C .....  INCREMENT  PARAMETER 

PhISAv*PhI  ( I ) 

Pm  I ( T )«AMAX1 ( AMIN1 (hTGN( I } . PHISA V.SDFL ) *LOW( I ) ) 

C,,t, .EVALUATE  FUNCTION 

UPNl*u(toHE*PL»l *PhT ) 

N«N  ♦ 1 
N3*n3. 1 

C OUTCOME  an  improvement  f 

IF  (UpHl  .LT.UPmO)  GO  to*? 

IFCN.GT. LIMIT)  GO  To  1 

C.,,,, INCREMENT  PARAMETER  IN  OPPOSITE  DIRECTION 
S(I)«-S(T) 

PhI (I )*AMAX1 ( AMIM (hIGM(I) .PHISAV.SDEL) »LOw(I) ) 

C.,.#, EVALUATE  FUNCTION 

UPMI*u(6HEXPLR2»PHl ) 

NaN.  1 
N3*n3t l 

c .... .outcome  an  improvement  ( 

IF (uPh! .LT.UPhO)  GO  to  2 

C.,., .RESET  parameter 

1 PnI (I )*PhI3AV 

GO  TO  3 

c.... .retain  8FST  FUNCTION  VALUE 

2 upmO*uPhI 

3 IF(N,GT,lImIT)  return 

« CONTINUE 

RETu»N 

END 

SjaHOuTINf  PATmv (UPhO.PhI ,PhO) 

DIMENSION  PMTdJ.PMnd) 

common  /B/R/  k .DElTa.ALPHA.l IMIT.N»N2»N3»EPSLDN,ThFTA(15) t 

* 5(15)  *PFPS»  trace  •LOW  f 15)  » HIGH  (15)  .PANGEd  5) 
LOGICAL  TRACE 

REAL  l 0* 

M*1 

C,,,,. obtain  Exploratory  increment 
DELTAaSQRT (OElTa/K) 

c .... .EXPLORATORY  INCREMENT  too  Small  ( 

1 IF (DELTA. lT.EPSlON)  RETURN 

C ....  .Evaluate  function 

2 UPMl*u(5MPATMV,PMl) 

N a N ♦ 1 

1 

TF(N,GT.l  IMIT)  RFUtpN 

C CONDUCT  EXPLORATION- 

CALL  E*PLR(UPhI .Ph! ) 

C Outcome  an  improvement  r 

IF (uPmT ,GE ,UPhO)  GO  TO  U 

z parameter  change  significant  ( 

DO  3 I * ) » i* 

7F  ( ARSfPwI  ( n-PMO(  I ) ) ,GT  , 1 .F-6TARs(p*0(  I ) ) ) GU  TO  10 

3 CONTINUE 
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u TE(N.GT.LIMin  RETURN 

GO  TP  (5*7*9)tM 

C.,,.. REDUCE  EXPLORATORY  i NCRE  ME  NT 

5 D£  L T A a 0 , 5*OE  L T A 

C DEFINE  PROJECTED  POINT  NFARpR  «a3e  POINT 

on  6 i * l • k 

t>  Pm  I ( n*HOUND(  I *PHO(  I )f,5*THFTA(I)  ) 

GO  Tn  \ 

7 “« 5 

C...,. DEFINE  PROJECTED  POINT  IN  OpPOSTTp  DIRECTION 
00  8 1*1  »K 
S(I)*-S(I) 

8 PMT(naB0UNDCI*PM0(I)-,5*TMfTA(n) 

GO  TO  ? 

R RETURN 

c..,, .RETAIN  best  function  value 

10  UPMO«uPhI 

RETURN 
ENO 

• • 
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norm 


This  S‘I8-PI»0GRAM  IS  CALlFO  P734 

IT  IS  USED  aITH  TESTS  TO  OPTIMA  *ITH 

..»..PwvslCAuLY  REALIZABLE  niELEcTRTC  MATERIALS. 

REFER  70  SUBROUTINE  EXPlR  FoR  ThE  TABLF  SEARCH, 

Subroutine  razor (Pwo»ep3min.rhO,eta»kappa) 

DIMENSION  Phi ( 15) »PhOC 15) 

COMMON  /RZR/  K,DELTAtALPHA.LlMlT.N.N2tN5.EPSL0N,THFTAC15) . 
• SCI  5) .FFPS.TRAcE»LPHf 15! # HIGH C 15 ). RANGE ( 15) 

I OGIC AL  FINISH, TRACE 
B£  A l LOa 

CALL  SECOND (TIME! 

Na  1 
N 1 » l 
N2«0 
N3»0 
Nu«0 

ic«o 


C EVALUATE  FUNCTION  at  STARTING  POTnT 

UPH0KUC6HRAZ0R1 , PHO) 

IF(TRACE)  wR I TE ( b * 1 I UPHOf (PHO(I) »Ial iK) 

1 FoRMATClHOtUXibhuPHo  »»G?0,10»15X,5HPHO  ■ , 3G20 , 1 0/ (1 1 X » 5G20 , 1 0 ) ) 

c... ..finish  criteria  satisfied  f 
IFCFINIShC UPWO ) ) GO  TO  14 

no  2 i * i »k 

2 SCI)*1,0 

EPSlON"tPSMIN*(ETAMKAPPA) 

C CONDUCT  PATTERN  SEARCH 

IFCTRACfc)  ARTTFCbftoO) 

100  roRyAT(7H  PAT3R1) 

Cull  p a T SE  R ( PHO  * UPHo ) 


c..,, .FINISH  CRITERIA  SATISFIED  r 
IF IFINTSH(UPHO) ) GO  TO  14 
DO  10  J ■ 1 , KAPPA 
N 4 a N • N 4 

UPRO»uCbH«AZOR?fPHOl 

IFfTRACt)  rtRlTECb.jI  N4»UPH0,fPHp(I),I«l»K) 

3 F0RMATC4H  N a, 15, lbX,6HUPHD  a,G?0.10/6X,5HPHO  ■»3G20,10/ 

* C 1 1 X , 5G20 , 1 0 ) ) 

Mil  s»N 

IF CN.GT. LIMIT)  GO  T0  1? 

C OBTAIN  NFa  STARTING  POINT  AnD  ExPl ORATORY  PARAMETERS 

delt  AaO ,0 

no  4 I a 1 , K 

RANDOmb«!,0+2.0*RANfC0.0) 

PHl(I)sBOUNDCI»PHOfl)tRHO*BANDnM*EPSLON/lOO.O) 

Tt CPHI f II -PhD c I ) ) /Range  Cl) *100.0 

4 DEL TAsDELTA*r*T 

DEL  T AaSGRT (DELTA /K  ) 

evaluate  function  amo  conduct  pattern  search 

UP*I«UfbHRAZnP3»PMT) 

N a N ♦ 1 
N 1 a \ 1 t 1 

IFCTBACEl  h»lTE'a»Sl  i*PH  J . fPHl (!) ,1*1 ,K) 

5 FORMAT ( 1M0*4X, NHUPHy  a , G? 0 , j 0 , 1 5 X , 5HPH I a,3G20.l0/ 

* Cl! X .SG?0 .10)) 

EPSlOnbEPSI  DM/ETA 
Tf(TPACb)  ^RTTE(OflOl) 

101  FORMAT(7h  PATS»2) 

Call  PATst°(PHl , *ph i ) 

C OuTCOhF  AN  IMPROVEMENT  [ 


C-2 1 


T F CuPh! .lT.UPHO)  GO  TO  7 
U <n,GT.LImIT1  GO  To  12 
C.,,.. PROJECT  N Ew  POINT 
L 1 A aO  . 

00  b 1*1  »* 

TbPhO(T )-PHI Cl ) 

T 1 cT /WANG!  ( I 1 *1 0 1,0 
TuFTA ( J )«T 

PHTCI)«BnUN0CI»P^O(l)»T) 

6 OtlUrnfcLT/WTltTl 
0 P m I a i j P H 0 

C conduct  p A 7 T F WN  ^0Vf 

iFfTWACtl  I TE  (6»  102) 

102  fodmuj (7h  p a TM V l 1 

CALL  R A T m y (c'PhI  »Ph7  ibhO) 

C OiiTcOmF  an  IMPROVEMENT  [ 

TpCUPhI.GF.UPHO)  GO  TO  9 
C WfTAlN  BEST  POINT  4N0  FUNCTyON  VA^-lfc 

7 PPmOSliPhT 

C PR0JFCT  NEW  POINT 

OF  L TAsO. 
no  B i a 1 • k 
TsPhI Cl)-PHO(J) 

T i *T  /RANGE CI)*10«.0 
ThFTA(! ) z T 

0-n ( I :*ppi ( n 

PhT CI )*BOuNDfI #PM0(I )♦!) 

0 OElTAsOEL TA+T1»T1 

• F fN.GT.L I ^ I T ) GO  TO  12 

C CONDUCT  PATTFRN  M0Vf 

IF  (TRACE!  oPITF(fa*lo5) 

105  F n » M A T ( 7 W PATMV21 

Call  pa tmv (uphI iPhi , Pho) 

C OUTCOME  AN  I MpRO VE  N T ( 

TF  ruPnl  .1 T ,UPH0)  GO  TC  7 

C FINISH  CRITERIA  SaTj SPIED  r 

B IfCF TNTShCUPHO))  GO  TO  ju 

10  CONTINUE 

CAUL  SPCONO(TlMtl) 
r!yEsT  rvF 1 • T T Mf 
*piTEcfe.i n n,tive 

11  FpRMAT(*ONO  CONCLUSION  6V  PaZOP  Af  TER* . I U i * FUNCTION  EVALUATIONS*. 

* * AN0*,F7,?,*SFC0nDS*) 

Go  TO  16 

12  CALL  SECONO ( T 1 m£ 1 J 
T IMt a T T 1 • T TMf 

VA. p 7 T e ( 6 . 1 51  T I me 

13  FnPMAT  CbRHOLTMI  T OM  FUNCTION  E V a t (JETTONS  EXCFtOED  R Y RAZOR  AFTER* 

* E7.?,RM  sECOnOSI 

GO  TO  ! 6 

1U  r AI.L  SFCONOCTIMt  1 ) 

T I ME  3 T I MF 1-TlMp 
wBTTFfO.ISl  Time, N 

lb  F0PvATC33h0CONVE»GFnCE  OP  T A y NE  0 By  HaZOR  IN.F7.2, 

* 1UM  SECONDS  AFTEP,yS,?1  H FUNCTION  EVALUATIONS! 

16  N u x N • N £l 

y R TT  E (6,171  N 1 « N2 i N3 , nu . UPHD • (phO(I).Iei*K) 

17  FORMAT  ( 3 7 h The  FUNCTION  F V A [ UA  T y Om  BPF  A K OO^N  IS/IlO, 

* Bh  By  RAZOM/I  1 0 . 1 Bh  p Y PfilMV,  AND/110.BH  BY  F*P!.R/ 

* 2«hOFTNAL  PAPAMFIfR  vAtUp.S.RH  N 3.J5.3*, 
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O n 


• oHt?HPi  s , G 1 5 . 7 » 5X  , S^PHO  • •r,15.7/(lSX*6G15.7)  ) 

R£  TORN 

end 

l C G X C A L FUNCTION  FInTSH(UPHo) 

COMMON  /RZ«/  K, DELTA. ALPHA. I I M I T , m . N? , N3 , E PSL ON , THE T A ( 1 5 ) , 

* S(l5)  ,FEPS.TRACE,L0H(15) »HIGH(15) ,RANGE(15) 
LOGICAL  TWACF 

BPAl  10* 

F I N I Sh*  , F AL  SE , 

IF  (UPhO.LT  .FFPS)  F I N T S«*  , T»l)E  . 

IF(TBACE)  *RTTE(6,100) 

100  FPRmaT(7h  FINISH) 

BpTuPN 

e nd 

subroutine  patser (Pmo.upho) 

DIMENSION  Phi ( 15) ,PHOfl ) 

COMMON  /RZR/  k, DELTa,6LPMA. LIMIT, N, N2»N3»£PSLON, THETA(lb)  , 

• 3(15)  .FFPS,TRAcF,LOWf 15)»MIGH(15),PANGE(15) 

logical  FINISH, trace 

RFAL  t0« 

1 I F ( N ,GT,  LIMIT)  RETURN 

C , , , , .exploratory  increment  too  small  ( 

IF(L)ElTA.lT.EPSLON)  RETURN 
IF (E INISH(UPHO) ) RETURN 
UP* I SUPHO 

on  2 I ■ i , k 

2 PhI  (I )«PmO(I ) 

r C ,,,,  .CONDUCT  EXPLORATION  AROUND  RASE  POINT 

* C *.t.L  fcXPLR  (UPHI  ,PhT  ) 

IFCTRACE)  HRTTE(6,100) 

ICO  FORMATf7H  PATSFR) 

C .... .OUTCOME  an  improvement  ( 

IE(UPhT.L T.UPhO)  GO  TO  3 

REDUCE  EXPLORATORY  INCREMENT 

7 DEL TA«ALPHA*OELT A 

GO  TO  1 

RETAIN  best  POINT  And  FUNCTION  VALUE 
project  NE h point 
DELI A*0,0 

'JP*0*UPH  I 
DO  4 I * 1 , K 
T 1 3 P H I ( I ).PnOf  I ) 

T»T 1 /WANGt  ( I)  * 1 00 , 0 

THE) A ( T ) * T 1 , 

pmp ( i j =p*i  ( : ) 

°h T ( I ) «B0UND ( I .PHOf  I ) tT  1 ) 

U DELTA*DtLTA,T*T 

IF  ( E I NT  SH  ( LlPHO ) ) RETURN 
IE(N,GT.L  I m I T ) RETURN 
C , . , , .CONOUC T pattern  mqve 

Call  patmv (uphi .phi ,php) 

TFfTPACE)  wRTTFtb.tOO) 
c ,,,, .outcome  an  improvement  t 

TF  (uPHl-HPH0)3,'»,7 
END 

FUNCTION  B0UND(  ! .fcxpR) 

COMMON  /p/R/  f.dELTa.ai  PHA  *l  IMT  t ,A',N?,N3,tPSLON.  Theta  (1  5)  » 

* S(15).FEPS, TRACE. L OH M5)»HIGH( 15) .RANGE (15) 

RfAL  l*T^ 

ROUNDa  am  A XI  ( AMTN1  (m  jGm  (I  ) ,F  *PR  ) ,1  0*(  T ) ) 
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Of  Tu»M 

FND 

Sl'HKOuT  INf  fc  XPLN(uPhD,PhI ) 
dimension  pmici) 

CO**PN  /» ZR/  K .DE  lTa  , ALPHA  T^TT  .K,»N?tN3.fc^3L  on,  THET4  ( 15)  # 

• 5(15)  »FEPSiTRAcF.lP*r  1E>)  *HIGN(15)  .pANGEUS) 

LOGIC  At  TRACE 
REAL  LPh 

COMMON  /I'E  C T / 0(lJ)tFPS(la),MSfASl  AflS.IP 
DIMENSION  Tarl(92) 

0 a T 4 T 6 HL.  /I .02. I .05, 1 .05, 1 . 1 . 1 . la. 1 . 15, 1 ,2. 1 ,23» 

• l,J*l.«*l.«l*l,5,l,6,llT,i.B,l,9, 

* 2,0,2. 1 «2,2,2,3«2,32»2.a»2.a5,2,a8»2,5»2»5?. 

* 2. 53. 2. 5a, 2, 55. 2. 62*2. 7, 2.8. 

* 3. D, 3. 32. 3. «. 3. S» 3. bb, 5. 75. 5. 78*3.78, 3. 82, 3. 8a, 

* U , 0 * a , 3 . u , a . a . 5 , u , 52 i 5 , 0 , 5 . 2 • 5 • B » 5 , 5b  * 5 , 65 . 5 , 68 , 

♦ 5,8,b,0»b,07,b,l,b,3.b.5.8,b*b,79,b,8»b,9, 

» 7,0,7.12,  8. 0.8. 2. 8. 5* 8.79,®.«»  9 . 0 , 9 . 5 , 9 . b , 

• 9.7, 9. 8,  10.0,10.1,  1 1 .0,  1 2 , 0 , 11. 0.  ia,o. 

* 15.0,  lb.O.  17.0.  18.0,  Iq.O,  20,0,  2l,0. 

• 22,0.  23,0  . 2a. 0,  25.0  / 

DATA  NB/q?/ 

DO  B 1*1.8 

5DFL*0ELTA/100.0*S(I l*RANGf (!) 

C INCREMENT  parameter 

Pw  T S A v*Ph  HI) 

PmT ( I ) BAM A XI (AMIM (kTGHf I ) .PHISA V^SDEL) ,L0W(I ) ) 

C.,,,,FOR  I . Lf , IP  »PHI  COMPONENTS  apf  0IFLECTRIC3* 

C ?cA«Ch  TARL  FOR  TmE  CLOSEST  niFtFrTklC  CONSTANT 

C.,,,.AND  REASSIGN  That  VaLUF  INSTEAD  Of  TmE  COMPUTED  VALUE 
T F f I , GT , IP)  GO  TO  9 0 
P^I 1*0,01 *AjNT (100»0*PHl (I  1,0.5) 
rsTvlNsqoq.O 
DC  88  L * 1 • N R 
TaBLLsTARL (l) 

T F ( hh  i I , E (3 , T A0L  l ) GO  TO  9o 
r'lST*A8S(PHl  I-T  A8LL  ) 

IPtDIST  ,GE.  D9TMIN5  GO  TO  88 
OSTmInbOTST 
L DC  *L 

88  CONTINUE 

Ph I(I)«TABL(LOC) 

C evaluate  function 

90  UPhi*u(6hE XPL R1 ,PHT ) 

NeNtl 
M J*N3* 1 

C . . . . , OijTCOmF  an  IMPROVEMENT  [ 

IF  t uPh! ,LT ,UPHO)  GO  TO  2 
TF  fN.GT.L IMIT)  GO  To  1 

C INCREMENT  Parameter  in  opposite  DIRECTION 

S(I)a-SU) 

PhI ( I )*ama* i ( AMIN! (hTGh( I ) tPHISAV.SOFl ) »lG"( ! ) ) 

FqR  T ,lE.  TP  *Pm I COMPONENTS  ARF  DIFLFCTklCS* 

SEamCh  T A HL  FOR  TnE  CLOSEST  DIFlFcTRTC  CONSTANT 

and  reassign  that  vaiuf  Instead  pf  tme  computed  vai uE 

IF (I  .GT,  TP)  GO  TO  100 

P“ I ; *G . v t « A 1 N T M oo . 0*PH  I f n ♦ 0 . 5 ) 

DStNTN*9PR.O 

DO  90  LSI »NR 

TA8LI sTABL (L ) 
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JF(P"II  ,EQ,  TA&Ll)  GO  TO  ! 0 0 
Dl8TaAH3f PHI I-TABLL ) 

IFCDIST  ,GE.  OSTEIN)  GO  TO  99 

OSTMI\aDIST 

10C*L 

99  rOMTIMIfc 

PHI ( I JsTABL (COC) 

C • • . t .F  V AU'ATE  FUNCTION 

100  HPhI«u(6HEXPLR?*PHT  ) 

N jN*  1 

N5»N3*  1 

C • . . . .OuTCO^F  ft N IMPROVEMENT  ( 

TFOjPhI  .lT.UPmO)  GO  TO  2 

C RESET  PARAMETER 

1 Pwl  ( I ) bPhISAV 

GO  TO  3 

.....RFTAlN  bEST  f uNC  T I On  value 

UPHO»uPhI 

IF  (N.C.T,  L!mIT)  return 

CONTINUE 

RETURN 

END 

SUBROUTINE  patmv (uPhO.PhT iPmO) 

DIMENSION  PhHU  .Pho(1  ) 

COMMON  /p z»/  K.DELTa*AI PHA. LIMIT. N.N?tNJ.EP8L0N. THETA(t51 t 
• S C 1 5 ) »FEP8. TRACE* LOW ( 151 »MIGH( 15) .RANGE ( IS) 

LOGICAL  trace 
REAL  lOi" 

M*1 

OBTAIN  EXPLORATORY  increment 

OELTAbSQRT i'DELTA/K) 

Exploratory  increment  too  small  r 

TFfPfLTA.LT.FPSLOM  RETURN 

Evaluate  function 

UPm i su  f 5wp aTmv  * PhI  ) 

N«N.  1 
N?*N?  t J 

• IFfN.GT,  LIMIT)  RETURN 
c .... .CONDUCT  FxPlORATTON 
CALL  EXPIR(JPhI.PHI) 

C .... .outcome  an  improvement  i 
IF (uPhT .Gfc.UPHO)  GO  to  u 

c..., .parameter  change  significant  i 
on  3 1*1  .k 

TF(ABS(PhI  n ).PHO(T  ) ) ,GT,1  ,f-6*ABSCPmO(D)  ) GU  TO  10 
3 CONTINUE 

« TEfN.GT. LIMIT)  RETURN 

GO  TO  (5.7*9). M 

C REDUCE  Exploratory  increment 

5 OFLTA*0,'ft*Ofl TA 
M*<» 

c DEFINE  PPOJECTEO  POINT  NFARFB  RaSf  POINT 

DO  b I * 1 . * 

6 Ph J ( I ) a HO  iND ( I « U01J).,5*THFTA(I)) 

GO  TO  1 

7 m,  j 

C.....DEF  ! N1  -»()■,  ' '■  •'■Ir  N OPPOSITE  DIRECTION 

on  « I a i • * 

S ( T ) s»S ( I ) 

8 PhI ( I ) shOunq  f I . phot  I ) -,5*ThfT a ( p ) 
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50  TO  ? 

« f»£Tu»N 

c bfst  EuNC  t i 0 n VALUE 

10  UPhO«uPhI 

PEtUPN 
E\n 

09 


f 
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prOGRa”  E I L TF w3 ( INPijT  .OUTPUT.TAPF  i » I NPUT  »TAPE6>0UTPUT  » T AP£R ) 
C..,.,ThTS  IS  The  OPI&INA(_  working  PLOTTING  PROGRAM 

01 ^ENS I ON  T(i 03) *ST( 1 03) .0(21 ) »EPS<22) *DIELCT(6) .FREQC6) . 

* SAV£(6*10l) 

COMRLEX  Z 

LOGICAL  ODD ♦ 8K I P ! T 

ntT a ST  /100*0,0  * ,999  . 0.0  . ,?0  / 

DATA  T /101*0.0  . .50.0  . JO.O  / 

HaTa  T^OPI  /fe,283l«530T179*/ 

CALL  PLOTEf30t«) 

00  100  1*1*100 
100  8T(l)*0,01*(I-i) 

INPUT  ns.  t hF  Number  of  NON. AIR  DIELECTRICS 
N SL A R S IS  ThE  TOTAL  NUMBER  OF  Ajr  % NON-AIR  OIElFCTRICS 
IP  IS  The  NUMBER  of  DIFFERENT  NON-AIR  dielectrics 
1 10  PRINT  *, I • 

print  • nr  slabs* • , 

PE  AD  ( 1 • * ) NS 

IF (NS  ,LE.  0)  GO  TO  190 

N Si ABS*2*NS 

N3L API »NSLA0S-1 

I P * ( N S ♦ 1 )/S 

C TnPjT  OIELCT(I)  1*1  To  IP  X INITIALIZE  EPS  VECTOR 

PRINT  *,  > E PS ( 1 TO  • ,IP.» )e«  , 

RE 4 0 ( 1 * * ) (OIELCT ( I ) , 1*1 , IP) 

ODD* . TRuF , 

if(mod(ns*2)  ,fq,  0)  odd*. false. 

00  ze  1*1  . IP 
Tr->2*I 

f PS ( I X • 1 )*1  .0 
E PS ( I X ) *0  I ELC  T ( I ) 

IF (ODD  .AND,  (I.EO.IP))  GO  TO  7p 

T x*2* ( TP*  I 1 

EPS  ( I X - 1 )»1  ,0 

I Y *NS / 2 * 1 -T 

EPS(  Ix)*OIFLCT( IV) 

76  CONTINUE 

C.....INPUT  0(1)  1*1  TO  f-SLA»l 

PRIM  *,'D(1  TO  ' *N8LA«1  , • )i'  ♦ 

»£  A o ( 1 . * 1 (Of  I) *1*1 .NSLAPl) 

C,,,, .INPUT  T f or  T m MODE  ....  * 

112  PRINT  *,fMODF*TE/TM*et, 

RE  AO  C 1 * 390 ) I NMQOE 

J90  eoBvat(a?) 
mqOE  *2wTm 

T F ( I NmODF  ,EO,  MODE)  GO  ’0  il3 

MOOfcs^HTF 

IF ( I NmOOF  ,EO,  RQDF)  GO  TO  1 1 3 
GO  TO  11? 

C.,,,  .INPUT  SIN(TmfTAI)  SIN(THETA2)  Dpi  Dfa? 

113  Print  ♦. »SThi »5lH2,nRi ,DB2#f , 

RE  A D ( 1 . * 1 SThj .STH2.0B1 .0B2 

r INPUT  FRFG  lOw. hIGh, incr 

PRINT  IO.mjGh,  INCR*  ' . 

PF  Al  ( 1 . * ) F l 0 » E h ] * F I N C P 

NFH«l!.T(  (FHI-FLO)/FTNfR  ♦ 0.5)  ♦ 1 

c riMtPM\r  if  data  printing  desiwfd 

litt  PR  I M »,I<ANT  PRINTING  DF  D4TA  •V/N*X<* 

Pf  AD ( 1 • 390  ) INTwri 
SkIPITs.TR.JE. 
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TFCINTRC?  ,tO.  1*N)  GO  TO  1 1 5 
SKlPITa.f ALSF  , 

TFdMTRCE  .EC,  1*Y)  GO  TO  1 l *S 
GO  TO  11 o 

US  tail  A X I s C t ,5.3,5.  low  SINO**i*-.-l  O.S,  *0.  *0.  . ,?) 

CALL  A x 1 s C 1 ,S.  3.5.?RHTRANSmtSSTON  COEFFICIENT  ( DB  ) . ?R . 5 . * 90  . , -50  . 
* .10,0) 

CALL  PLOin  .5, 3. 5.-3) 

L s 1 

DO  *30  I a 1 » NF P 
FbFlO+  C I - 1 ) AFINCO 
TPIL=TwQPI*F 

FREQ (L ) «F 
L»l  ♦ 1 

DO  120  Jel * 1 0 1 

CALL  R*TX(NSlA6S.ST(J)fTPlL,*0DF.FPS.D.Z) 

TJ*CA8S(7)**2  t 

tJb10,*ALDG10(TJ) 

T (J)»AMAX1  f TJ.-50.0) 

SAVE  t I . J ) «T ( J ) 

120  CONTINUE 

CALL  LINE  (ST  ,T  . 1 01  ♦ i .S.  I»P 
130  CONTINUE 

C TRACE  OUT  SAVED  n A T A 

aPTTF (b,3«U)  (SAVE  (K  » n »K«1  ,NFP) 

3R-J  FPRvi T ( 1 X .0 (F7, 3. 3X  ) ) 

IF (SKJPIT)  GO  TO  30c 

00  200  J » 1 . 1 0 1 

•'  < I T F ( 6 . 3 R 3 ) S T ( J ) . f S A VF  ( K , j ),  K s 1 , NF  R ) 

3R3  FORMAT f lX.F5,J.2X.b(F7,3»3X) ) 

200  CONTINUE 
JOC  S = 5 , *STh1 

CALL  PLOT ( S I 0 , o » 3) 

CALL  PI OTCS.5,5.2) 

Call  SYMBOL (R9R. .5,b.,07.aHgTHi,o..U) 

SbS,*STH? 

CALL  PLOT ( S i 0 , 0 ,3 ) 

CALL  Pi OT (S.5,5.2) 

CALL  STMBOL(RRP..b.6».07.aHSTH?,O..U) 

SaS.*0. 1 *081 
Call  PLOT ( 0 , # S « 3 ) 

Call  plCTCS, .8.2) 

CALL  SYMBOL  (5.  1 ,RR9.  . . 07, 3*D«1 .0.  ,3) 

3*5, *0. 1*0«2 
CALL  PL0T(0,,8,3) 

CalL  PLOTCS.  .8.2) 

Call  Symbol ( 5 , 1 . .07 • 3*082.0, ,3) 

Call  SYMBOL (1 , ••1.5..1 .6HLFGFND.0, .6) 

1 *0 

Y * • 1 ,5 

no  1 3S  1*1.5 

Y • Y ■ , B 

no  135  J*  1 .2 
l *1  * 1 

X«.5*2.5* ( J*1 ) 

TF(L  . r,  T , NFR)  GO  Tn  13b 
CAlv.  symrolcx,v,.i»l-i  .0..-1) 

Call  5 ympOL ( X * ,2 . Y . . 1 ♦ THFHFo  b ,0..T) 

CALL  Number  ORB. . RQQ. »,1«FRf0(L)«0.»?) 

135  CONTINUE 
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136  CALL  PLOT ( T , , » 3 , S * •%} 

GO  TO  110 

UO  CALL  PLOT:0.,0,*R9R) 

endfIlE  « 

STOP 

f* 

V* 

END 

SuBROuT I M E PXTx(NSLaR«SKN*TpIl*mOd£«FPS*D»PT ) 

COMPLEX  PT  * U * J * GLR  T ,TY*S(2*25  • A ( 2 , 2 ) * AD ( 2 1 *2*2) 

DIMENSION'  FP9Cl)tD(iJ»Y0(21).Aj(2i»2,2) 

0 E A l KAPPAC21) 

EQUIVALENCE  (SU.Sd.in  . (821.9(2*1))  . (PLTY.TY) 

DATA  J/(0.  • l . )/ 

C 

C COMPUTE  KAPPAVS  AND  Y 0 V S ARp  AL* 

c 

9kN2b3kn*Sxn 

DO  100  NDbI.NSLAB 

KAPPA  C NO)* 3QPTCEPS(ND)-SXN2j 

VO(ND)skaPPA(NDJ 

I F ( MODE  ,EO,  2HT*)  yO ( NO ) .EPS ( Np ) /KAPPA ( NO ) 

100  continue 
C 

C . , , , .CONSTRUCT  JUNCTION  wAVF  TRANSMISSION  MATRICES*  Aj,  aREAl* 

C 

N8lA81«NSlAB»1 

DO  no  NDe  1 * N3L  ABi 

PLTY»l.O/(YOCNO)6YO(ND6lJ) 

911sfY0(ND)-Y0(ND*l)l*RLTY 
32;s2,0*YO(NO)*PLTY 
5LtTcY0  C NO ♦ 1 ) / Y 0 ( NO ) 

AJCNO* 1 »l)«i .0/321 
AJfND, 1 *21*31 1/321 
A J ' NO , 2 . 1 ) *A J (NO  * 1 ♦ 2 1 

AJfN0,?«2)B(RLTY*S2l*S2lYSU*8in/92l 

C 

C. ...  .CONSTRUCT  LINE  length  wave  TRANSMISSION  MATRICFS.  AO.  AC OMPi_F  X * 
C 

AOfND|1*l)»CFXP(J**<APPA(NO*l)*TPTl*D(N01) 

A0CNn,2*?i«i./A0(N0,i«n 

AD(ND* 1 .2)8(0,  *0. ) 

110  ADC N 0*2*1 ) » ( 0 , . 0 • ) 

L 

C FOR*  total  WAVE  ThanS^ISSiOn  matRtX*  A,  ACOMPlEX* 

C 

DO  1 3 0 NDr 1 ,? 

DO  130  MDc 1 ,2 
A (ND«MD)a(0. «0„  ) 

IF  (NO. EQ. MO)  ACND.mO)»{1,.0.) 

130  CONTINUE 

c 

DO  ISO  NDal *NSL ABl 
NOt  *NSl AR-ND 
DC  160  I Da  1 *2 
OH  160  K Da  1 .? 

Ty*(0, .0.) 
no  150  L0>*  1 .2 
U*f0.,0.) 
no  i**o  mds  i ,? 

l«0  U«U*AO(NDl . LD»Mo3 fA (My.wO) 
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SYNTHESIS  OF  PLANE  STRATIFIED  DIELECTRIC  SLAB  SPATIAL  FILTERS  U— ETC(U) 
DEC  76  J H POZGAY#  S ZAMOSCIANYK*  L R LEWIS  F19628-76-C-0189 


BR-9389 


RADC-TR-76-408 


00K: 

|kt — 

IF" 

Jt^'l 

*r 

• 

DATE 

FILMED 

•5-77 


I 

150  TY*TY*U*AJ ( NO  1 * 1 0 ♦ L D ) 

160  3(I0»KD)»TV 

00  170  I0»!»? 

00  170  KO«lt? 

UO  A(ID*K0)*8(I0tKD) 

180  CONTINUE 

GLRTs(V0(N3LAB)-Y0(l ))/CVO(NSLAR)fVO(l )) 
TYbA(1«1)*G|.PT*A(1«?) 

C RETURNS  TRANSMISSION  COEFFICIENT  ^COMPLEX* 

RT*( 1 .♦GL«T)/TY 

RETURN 

END 

00 
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oor*  •—  niSnnnn 


PROGRAM  F I LTR4( INPUT* OUTPUT, TAPfl .INPUT »TAPC6«0UTPUTf TAPER) 

MODIFIED  FILTERS  PROGRAM  TO  WORK  0N  6 SLAB  CASE 

.....WITH  SET  STHl .STM2,FL0.FHI.FI^CP.DH1.0B2 
.....ONLY  DIELECTRIC  ♦ THICKNESSES  ARE  INPUT 

these  input  values  are  plotted  on  the  graph, 

.....ROTE*  ALL  AIR  DIELECTRICS  HAVE  BpFN  ASSIGNED  EP3*1,02 
,,,,,  TO  CONFORM  WjTH  A PHYSICAL  FOAM  APPROXIMATION  TO  AIR, 

DIMENSION  T ( 1 03 J .ST(i03).D(2l)*EPsf22)*OItLCT(6).FRCGC6). 

* SAVEC6.10D 

COMPLEX  Z 
LOGICAL  ODD 

DATA  ST  /100W0.0  .*  .999  , 0,0  . ,?0  / 

DATA  T / 1 0 1 *0 . 0 . .50,0  , 10.0  / 

DATA  T*OPI  /6. 2831853071796/ 

CALL  PLOTF(30.«) 

DO  100  1*1.100 
00  ST(I)*0.01*(I«1) 

8THi«, 17a 
STH2«,U2 
DR  1 *•  1 • 0 
DB2»-10,0 
FLO*. 98 
FhI*1 ,02 

fincr«.o? 

NER*INT(  (FHI«FLO)/eINCR  ♦ 0.5)  ♦ 1 
modE*2HTm 

.....INPUT  NS,  THE  NUMBER  OF  NON. AIR  DIELECTRICS 

NSLABS  IS  THE  TOTAL  NUMBER  OF  A I R X NON. AIR  DIELECTRICS 

IP  IS  the  number  OF  DIFFERENT  NPN.AIR  DIELECTRICS 

NS*6 

N3LA83«2*nS 


nslabi*nslabs-i 
IP*(NSFl )/2 

c..., .INPUT  DIELCT « I ) 1*1  TO  IP  X INITIALIZE  EPS  VECTOR 
110  PRINT  * , * EPS ( 1 TO  ',IP.*)«», 


ad c i . » ) (dtelct(I).i«i.ip) 

IF (DTELCT ( 1 ) ,LT.  0.0)  GO  To  i«0 
ODD*, TRUE. 

lE(MODfNS.2)  .EG.  0)  ODD*,FaLSF. 

DO  78  1*1. IP 
lx«2M 

EPSCIX-l )*1 ,0 
EPS ( I X) *DIELCT  fl ) 

TEfODO  .AND,  (I.EU.lP))  GO  TO  78 

I x*2* ( I P ♦ I ) 

eps(  ix-n*i  .o 

I y«NS/2t 1 - I 

EP8(IX)*0IELCT(IY) 

76  CONTINUE 

c input  D(i)  1*1  to  n9la»i 

PRINT  F , • D ( 1 TD  • .NsLABl , • )■' » 

Rf  AD ( 1 • * 1 (0(1) « !*1 .NSLAfll ) 

CALL  AXISC 1 .5.3.5, lQHSINF  ThET A 0 .5 • .0 . . 0 , . , 2 ) 

CALL  AXISM ,5»3.5*29HTRANSMjSSION  COEFFICIENT  (D8) » 29.5. . 90 J ,.50. 
* *10.0) 
call  PLOT ( 1 ,5. 3.5*. 3) 

l si 

DO  1 3o  IbI.NFR 

FsFlD. (1*1 )*E inc« 

TP1L*TwOPI*F 
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freq<l)»f 

DO  120  J«1 » 101 

CALL  RXTX(NSlABStST(J) ,TPlL,MODp,fPS,D#Z> 
TjaCABS(Z)**? 

Tj«10.*ALnG10(TJ) 

T ( J ) «AMAX 1 f TJ#-50.0) 

SAVE  C I ♦ J)aT(J) 

120  CONTINUE 

CALL  LlNE(STiTtl01.l»5»I-lV 
130  CONTINUE 

C TRACE  OUT  saved  data...  ...... 

WRITE(6»394)  (SAVE(K»l)tK«l.NFR) 

39a  fopmat(ix»6(F7, 3*3*1) 

S*5 , *STH 1 

CALL  PL0T(S.0.0,3) 

CALL  PLOT (S,5«5,2) 

CALL  SYMBOL (R99,«5,6« .07 ♦ aMjTMl , ft . , U ) 

S»5,*STH? 

CALL  PLOT ( S » 0 , 0 » 3 ) 

CALL  PLOT (S»5,5»2) 

CALL  STMROLC999,»5,6«,07taH8TH?.O.»ai 

S«5.y0.1*DB1 

CALL  PLOT(0,.S»3) 

CALL  PLOT (5. *S»2) 

CALL  SYMBOL (5, 1» 999., . 0 7 , JHDB l , 0 . , 3) 

S»5. ♦0,1*OB2 
CALL  PLOT(0.,3,3) 

CALL  PLCT(5,.S»2) 

CALL  SYMBOL (5,1 ,999., . 07 , 3H0B2, 0 . , 3) 

CALL  SYMROLd  ,*»'.75,.1,6HLEGEN0,0.#6) 

L«0 

Y»»0, 75 
DO  135  X«1*S 
YbV.,25 
DC  °1  35  J*1  »2 
L *L  ♦ 1 

Xa,5+2,5*(J*1 ) 

IP ( L .GT,  NFR)  GO  Tn  136 

CALL  SYMROL  (X,Y4.05,  .1  ,L»1  iO.fl  ) 

CALL  SYMBOL (Xf.2»Y,.l tTHPREO  • ,0.,7) 

CALL  NUMRE»(999.,990.,.l ,PRpO(L),ft.,2) 

135  CONTINUE 

136  CALL  SYM8OL(0,0,«1.75,,1,1?wDIElPcTPTCS«,0.,12) 

DO  200  Haiti 

CALL  SYMBOL (999,0 ,999.0, . 1 ,?H  ,ft.»2) 

200  CALL  NUMBp9(9 99, 0,999.0, . 1 » DIELC  T (K  ) ,0, ,2) 

CALL  Symbol (0.0, -2.0, . 1 ♦ 1 2MTHlCKNp3SFSa, 0 , » 12 > 

CALL  ftHERE (X0,Y0,F ACTR) 

C...  ..SINCE  amfaf  mas  BEEN  CALLED  AFTfR  A CALL  TO  SYMROl 
C.....TmE  COORDINATES  RETuBNFO  TO  X0»YO  AHF  THOSE  OF  THE 
C,.,,,LA3T  PART  DF  THE  SYMBOL  DRAmN, 

Co. ...hence.  They  REOUIRE  A SLIGHT  adjustment  to  LOCATF 
C.,.,.ThE  COORDINATES  OF  THE  BEGlNNJNr,  nF  THE  NEXT  SYMBOL, 
C.....TME  NEXT  LINE  ADJUSTS  FOB  ThF  • SYMROL  WHICH  *AS 
C..,,,ThE  last  drawn  BEFOrF  THE  CALL  TO  WHERE (XO,YO,FACTB) 
Y0*Y0-,05 
L«0 

DO  202  K s i , 6 
XaXO 
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Y»Y0-0.28A(k»1 ) 

DO  201  KK»t  ,4 
L • L ♦ t 

IF(L  ,GT,  NSLABl)  GO  TO  159 
CALL  SYMROKY.Y,  , V,2H  ,0..?) 

CALL  NUMRER(999,0»9q9,O. ,1 .D(L) .0..3) 

*■999,0 

Y«999,0 

201  CONTTNUE 

202  CONTINUE 

139  CALL  PLUT(7,,.3,5.*55 
GO  TO  no 

UO  CALL  PLOT(0,.0,i999) 

ENDFIlF.  a 
STOP 

c 

END 

SURPOuTlNf  HXTX(NSlA8»9*N,TPll.MOoF»EP3.D»BT) 

C 0MPLE X PT,U,J.GLBT,TY.3(2«2).A(2#2),AD(21.2.2) 

DIMENSION  EPS( 1) .D(l ) »Y0(21).AJ(2i .2.2) 

PfAL  KAPPAT21) 

EQUIVALENCE  (31 1 # 9( 1 * 1 ) ) » (821.3(2.1))  ♦ (RITY.TY) 

DATA  J/(0,.l,)/ 

C 

C.,.. , COMPDTF  KAPPAVS  ANO  Y0\S  PBf AL* 

C 

SkN2*SKN*sKN 

00  100  NObI.NSLAB 

KAPPA(ND)»SGPT(EPS(mO)«SXN2) 

Y0(ND)*KAPPa(ND) 

IF(wOdE  ,EQ.  2HTM)  Y0(ND)*EP8(Np)/KAPPA(ND) 

100  CONTINUE 

c 

c* , , , .CONSTRUCT  JUNCTION  WAVE  TRANSMISSION  MATBICES.  AJ,  PBfALP 
C 

N8LAB18NSLAB-1 

00  110  NO* 1 . NSL A8 1 

BlTy81,0/(Y0(ND)*YO(NO4D) 

81 t«r yO(N0)»VO(ND*| ))*PLTY 
821*2, 0*YO(NO)*RlTY 
PLTVSYO(NO-H  )/YO(ND) 

A J ( NO , l . ] )»l ,0/321 
AJ(n0.1 .?)*S1 1/821 
Aj(ND,?.l )*AJ(ND.1«2) 

AJ(N0.?.2)«(BLTy*82i*821YS1 1*S1 1 )/S2l 
C 

C.... .CONSTRUCT  UNE  LFNGTH  WAVE  TBANsmtSSTON  MATRICES.  AD.  PCOMPlFX* 

c 

AD (NO, 1 . 1 )*CEXP( J*xaPPA(ND^j  >PTPJi *D (ND) ) 

A0(N0,2,?)»1 ,/AD(NO, 1 .1 ) 

AD(N0»1 »2)*(0,»0,) 

110  Ap(ND,2.n*(0.»0,) 

c 

C.....FOPM  TOTAL  wave  TRANSMISSION  MATBtX.  a,  PCOMPLEXP 
C 

Do  1 3o  nosi,? 

on  i ?o  «n*i,? 

A(N().MO)«(0,  .0,  ) 

Tf  (NO, tQ. MR)  A ( Np , mD  )■(!,.{),) 

130  CONTINUE 
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c 


no  iso  no» i * nsl  ah  1 

NOl  SNSLAH-NO 
DO  160  I D * 1 t ? 

DO  160  KD«1.2 
TY*(0.»0.) 

DO  150  LD«lf? 

HafO.tO,  ) 

DO  1«0  mo«1,2 

1«0  UaUtADfNDl tLD*MD)*« (mD.KD) 

150  TY«TY*U*AJ(N01 *ID»LP) 

160  S(TD*kD)sTY 

DP  170  I0«1.2 
PO  170  KD« 1 » ? 

170  A(ID»KD)*S(ID.KD) 

160  CONTINUE 

CL»Ts(Y0(KSLAb)*Y0(n )/(YO(MSLAR)4YO(l )) 

T Y*  A ( 1 » ntCL»T*Afl  »?) 

C RETURNS  TRANSMISSION  COEFFICIENT  ^COMPLEX* 

RT«( 1 .♦GLRT)/TV 

RETURN 

END 


PROGRAM  CONTOR  (INPUT. 0UTPltT»TAPF5«lNPUTtTAPE6«0UTPUT) 
DIMENSION  EPS (22)  *D(21).U0(51).V0(51) 

100  READ  (5.800)  NSL AB8 

TF  (EOF(S).NE.O)  C*tL  EXIT 
NSLABaNBl  ABSM 

READ  (5.810)  (EPS(I). Isl.NSLAB) 

fP3(l>al. 

READ  (5.810)  (D(I)tial.N SLABS) 

DO  110  I ■ 1 .50 
U0(I)»f  I-1)*0,02*1.E-05 
V0(X)a(X«n*0.02Pl.Ea0S 
110  CONTINUE 

CALL  CROSS  (UO.VO. 50.50. NSLaB.EPS, D) 

GO  TO  100 

c 

800  FORMAT  (1615) 

810  FORMAT  (8F10.0) 

FND 

subroutine  cross  (uo.vo.nu.nv.nslab.eps.d) 

DIMENSION  U0(1).V0(1).EPS(1 ).D(1), INOM(51 .51). ICRS (51.51) 

COMPLEX  TTM.TVE 

DATA  TPIL/6.?83185508/ 

DO  110  IU»1.NU 
UaUO ( I U ) 

91 ■! •/( 1 • -U**2) 

DO  100  IVal.NV 
VaVO (JV) 

INDM(IU.IV)al00000 
ICRS(IU.IV)al00000 
IF  (U**2*vaF2.GE.1.0)  GO  TO  100 
WaSQRT(l.-U**2*VF*2) 

S2aSl/(l.-#**2) 

3TH«SuRT(l ,*w**2) 

CALL  RXTX  (NSLAB.STh.TPIL.2hTM, EPS* O.TTM) 

CALL  RXTX  (NSLAB»3Th,TPIL»2mTE»EP8*D.TTE) 

C 1 »0* V*w*S2 
C2eV*V*S2 
C3»S2*  (UP* ) **2 
T«CAB8<C?*TTM*C3*TTe> 

IF  (T.LT.l.E-10)  Tai,E-10 
To*200.*AL0G10(T) 

TbAMInI (T.999, ) 

TNOM(iU.IV)aT 
TaCABS(ClF(TTM»TTE) ) 

IF  (T.LT.l.E-10)  Taj.E-10 
T«.200.*ALOG10(T) 

T*AMIM  (T.99P.) 

ICRSdU.TVjaT 
100  CONTINUE 

no  continue 

write  (6.900) 

NI!JaMjN0(25.NU) 

DD  120  IVal.NV 

WRIT!  (6.910)  (INOM(lu.IV)tXllBl.NfU) 

120  CONTINUE 

MIDaMjNO  (51  » Nti ) 

WPTTE  (N.900) 
on  13c  IVsl.NV 

WRITE  ( 6 » 9 1 o ) ( JN0M(IU,IV)«lUa26*NlU) 
liO  CONTINUE 
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WRITE  (6.920) 

N I U»M  I NO ( 25 » Nu ) 

DO  1«0  IV«1.NV 

WRITE  (6.910)  (IC«S(TU.IV).IU*1,NT‘J) 
uo  continue 

wRTTE  (6.920) 

N JuaMjNO  (SI  »NU) 

00  150  IVal.NV 

WRITE  (6.910)  (1CRS(IU.IV)»IU»26»NIU) 

150  CONTINUE 

PfTOPN 

c 

900  FORMAT  ( 1 H 1 , 10X.7HNOMINAL  ♦/) 

910  FORMAT  (IK. 2614) 

920  FORMAT  (Ini  .1GK.5HCP0S9./) 

END 

SuBROljT  I Mf  RXTX(NSL AR.SKN.TpTl  .MOnE .FPS.D.PT ) 

COMPLEX  RT.U.J.GL«T,TY.$(2. 2), A(?,?),AP>  (21*2.2) 

Dimension  EPS(l).0(l)»Y0f21).AJ(?i .2.2) 

RfAL  KAPPA(21) 

EQUIVALENCE  (311*3(1.1))  . (321.3(2.1))  . (PITY, TV) 

DATA  J/(0,.l,)/ 

C 

C... ..COMPUTE  K A PP A \ S AND  YOVS  •PfAL* 

C 

SkN2»S*N*SKN 

on  10o  Nnal.MSLAR 

kappa (ND)sSQRT (EPS (ND)-SKM?) 

Yrt (ND ) BK AP®A (NP ) 

T F ( MODE  ,EO,  2HTM)  Y0(ND)aEP8(Nn)/KAPPA(ND) 

100  CONTINUE 
C 

c.,., .CONSTRUCT  tui^CTlON  WAVE  TR4nSmT$8tON  MATRICES.  AJ.  *RE  AL* 

C 

^SLA01«NSLA8»1 

DO  HO  NOsl.NSLABl 

BLTYal #O/(Y0(ND)*Y0(ND*l ) ) 

31  1*(Y0(MD)»Y0(NDM  ) )*PLTy 
321*2. 0*Y0(ND)*RLTY 
RL  T ya Y 0 ( NO* 1 ) / YO ( NO ) 

AjfND. 1 . 1 )*1 .0/321 
AJ(MO,  l .p)«3l t/321 
A J ( NO , 2 . 1 )»A.J(ND*1  *2) 

Aj(N0.?.?)a(RLTY*S2t*32l*S1 \*Sl  D/S21 
C 

C CONSTRUCT  LIME  LENGTH  wave  TRANSMjSSTOM  MATRICES.  AO,  6COMPLEXP 

C 

AO(N0.1 . 1 )*CFXP( J*KAPPA(NO*l )*TPTL*0(NP) ) 

AO(Nf). ?♦?)■!  ,/AD(ND.  1 . 1 ) 

AD(nO,1 *2)*(0..0.) 

110  AD(NO,?o  )*(0,.0.) 

c 

C..,,.FO»M  TOTAL  wave  TRANSMISSION  m A T R J X » A.  ^COMPLEX* 

c 

on  l So  NOal .2 
DO  1 It)  MOal  .2 
A (NO.MO)a(O,  .0, ) 

IF  (Nfj.En.MO)  A(ND.wO)a(1  ,,(»,) 

130  rOMTNOk 

c 
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DO  160  NQsltNSlARl 

NDl»N3L4fl-M0 

00  160  IOal.2 

On  160  KD»1 1 2 

Ty«(0#»0,J 

00  ISO  10>1»2 

l»a(0,«0.) 

00  160  MO«li? 

160  UaU*AO(NOl »LO»MD)*A(Mp#nO) 

ISO  TvarVtU6AJ(N01»10»LD) 

160  S(IO«KO)*TV 

00  170  IO«t  «2 
00  170  KDalt2 
170  A(IO«KO)aS(in«KD) 

110  CONTINUE 

6L®Ta(V0(iMSLAB)»Y0(n  )/(V0(iaSLA|)4V0(l)) 

T y»A ( 1 • 1 ) ♦6L*T*A ( 1 • 2 ) 

C...,.«EtURNS  TRANSMISSION  COEFFICIENT  aCOMPLEM* 
PT*C1.6CL»T)/TV 

Rf turn 
fno 
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